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Performance Report 


This project involved the installation of a downsized multichannel mass 
spectrometer instrument on the NASA P-3B aircraft and its subsequent use on the PEM- 
Tropics B mission. The new instrument performed well, measuring a number of 
difficult-to-measure compounds and providing much new photochemical and sulfur data 
as well as possibly uncovering a new nighttime DMSO source. The details of this effort 
are discussed below. 

The instrument used on PEM-Tropics B was the result of a previous downsizing 
project, and thus apart from some brief test flights the previous summer, it was 
participating in its first real mission. After the downsizer test flights several 
modifications were required prior to deployment on PEM-Tropics B These were largely 
completed in the fall of 1998 and some additional laboratory testing was also 
accomplished in this same time frame. The new instrument was shipped to the Wallops 
Island Facility in very early January, and was installed shortly thereafter, delayed 
somewhat by the lack of an aircraft. The instrument initially had to be configured to 
operate from only 60 Hz power and later switched to a combination of 60 Hz and 400 Hz, 
as aircraft power became available. Installation went smoothly, but ground-based testing, 
particularly of the DMSO channel, was hampered somewhat by the large amount of 
hydrocarbon background in the closed-up hanger. Test flights just before departure 
yielded good results for the OH, H 2 SO 4 , and MSA measurement channel but some 
problems in the DMSO/NH 3 channel. The flight from Wallops to Dryden showed 
improvement of the DMSO channel, but still not proper functioning. There was little 
time on the transits from Dryden to Christmas Island for repairs, and while the OH, 
H 2 SO 4 , and MSA channel continued to function well, problems persisted or worsened in 

the DMSO/NH 3 channel. 

The two inlets used to measure the above sets of compounds are very complex 
because of the sticky and/or reactive nature of the species. (For some of these 
compounds a simpler heated inlet might be used, except that such an inlet could 
evaporate the compound being measured from aerosols, thus causing a measurement 



interference.) Therefore, a physically complex but nearly wall-less ambient temperature 
inlet is required. In an attempt to work on components inside of one of these complex 
inlets (DMSO) at Christmas Island, a bolt was broken off, which precluded future access 
to the internal components including the ion source at Christmas Island. It was not until 
additional tools could be brought to meet the aircraft at Tahiti that the bolt could be 
drilled out and the inlet removed. With the inlet off, loose fittings on the ion source were 
found and repaired. In hindsight, a greater effort should have been made to gain access to 
the ion source at Christmas Island, but at the time it was not clear that the critical problem 
was inside the inlet, and tools and time were very limited at Christmas. Once repaired, 
the DMSO channel provided much new and exciting data during the flights out of Tahiti 
and on the return portion of the mission. The OH, H 2 SO 4 , and MSA portion of the 
instrument continued to function well throughout the mission, with the exception of one 
or two brief periods when flights through very wet clouds caused shorting problems in 
the ion source. At the end of the mission, the P-3 returned to Wallops, our instrument 
was unloaded, and was shipped back to our laboratory for further calibration and testing 
of both the OH, H 2 SO 4 and MSA and DMSO and HN 3 channels. The remainder of the 
first year was spent on data analysis. Overall, the new instrument performed well 
providing a great deal of new and interesting data. 

The second year’s effort on the project started just after the completion of the 
PEM-Tropics B mission. Data analysis continued in conjunction with more laboratory 
testing of the OH instrument. Improvements made to the down-sized OH instrument 
flown on PEM-Tropic B provided the first in situ calibrations of OH at high altitudes 
during daylight hours In the past, in situ OH calibrations were done only below about 
10,000 foot altitudes or at night when background OH was very low. The new high 
intensity dual output photolysis source used on PEM-Tropics B dramatically increased 
the number of calibrations made above 10,000 feet (because most flights take place 
during daylight hours when levels of ambient OH are as high or higher than those 
produced by the old calibration source at the low water concentration typically found 
above about 10,000 feet). These new calibrations lead to the surprising result that the OH 
instrument becomes more efficient and sensitive at higher altitudes. Initially not sure of 
this result, a laboratory study of several months duration was begun to verify and help 



understand these findings. By the end of the study, it was concluded that the initially 
observed results were correct, and that the increased sensitivity for detecting OH was due 
to improved efficiency for sampling OH at lower pressures and temperatures. The 
decrease in pressure and temperature at high altitudes leads to a reduction in the Reynolds 
number in the ion source region, which reduces turbulence, wall loss, and mixing and 
results in a higher percentage of ambient OH reaching the measurement region. Flow 
prior to the calibration source is not affected by these parameters because the central axial 
flow does not come in contact with the walls and the actual flow velocity is measured. 

The final OH H 2 SO 4 and MSA data were then submitted to the NASA data archive 

DMSO measurements were also conducted during PEM-Tropics B and resulted in 
some exciting new results. Most important among these is that DMSO appears to have 
an unexpected nighttime source. These measurements, however, also required many 
laboratory tests to insure their validity, since this was the first time that the ion chemistry 
scheme used to measure DMSO in PEM-Tropics B had been used on an aircraft. 
Calibration sources had to be retested and ground calibration had to be compared to in 
situ aircraft calibrations with isotopically labeled DMSO. The possible interference by 
several compounds had to be tested. The result of these tests suggests that the high 
observed DMSO is real, and no significant interferences could be found. The uncertainty 
in the data is, however, somewhat larger than anticipated due to problems with some 
commercially available permeation cell calibration sources of isotopically labeled 
DMSO. The DMSO data was then also submitted to the NASA data archive. 

Ammonia measurements were less successful. While ammonia can easily be 
observed by our instrument its very sticky nature prevented a meaningful measurement. 
Even with very little exposed wall area in our sampling and ion source regions of the 
instrument ammonia was still lost to the walls and later re-emitted so that the measure 
ammonia signal had little relation to ambient ammonia concentrations. 

The highlights of our PEM-Tropics B measurements were: an extensive mapping 
of the OH radical over the tropical Pacific Ocean, and the discovery that there must exist 
a nighttime source for DMSO (an important sulfur oxidation product). 

The results of the present work were presented in a number of talks at the special 
spring 2000 AGU session in Washington, DC. A list of these talks is included as 



Appendix A. This research effort also resulted in a number of manuscripts, which are 
listed in Appendix B. A more detailed discussion of the measurement techniques used 
and the scientific significance of the resulting measurements is included in the 4 papers 
included in Appendices C-F, which were written by our group. 
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Appendix A 


Talks resulting from PEM-Tropics B 

American Geophysical Union 2000 Spring Meeting, Washington, DC, May 30-June 3 

1 . Lee Mauldin . D. Tanner, C. Cantrell, E. Kosciuch, C. Gao, J. Nowak, D. Davis 
and F.L. Eisele. Hydroxyl radical measurements performed aboard the P-3 
aircraft during PEM-Tropics B. 

2. J. Nowak . D. Davis, G. Chen, F. Eisele, L. Mauldin, D. Tanner, C. Cantrell, E. 
Kosciuch, A. Bandy, D. Thomton, A. Clarke and V. Kapustin. A new source of 
marine sulfur? 

3. R. Weber, B. Anderson, L. Thornhill, D. Tanner, R. Mauldin, E. Kosciuch, C. 
Cantrell, F. Eisele, A. Clarke, V. Kapustin and K. Moore. Comparison of 
atmospheric physical and chemical conditions under which large marine new 
particle formation events occur. 

4. G. Grodzinskv, D. Davis, F. Eisele, L. Mauldin, D. Tanner, C. Cantrell, J. 
Crawford, S. Liu and Y. Wang. Latitudinal/altitudinal PEM-Tropics A/B OH 
distributions: Comparison of observations with models. 

5. G. Chen, D. Davis, L. Wang, J. Crawford, D. O’Sullivan, B. Heikes, F. Eisele, R. 
Mauldin, D. Tanner, C. Cantrell, E. Kosciuch, B. Anderson, S. Liu, R. Shetter and 
B. Lefer. Boundary layer HO x and O 3 chemistry based on PEM Tropics B 
airborne observations. 

6 . D. Davis. G. Chen, J. Nowak, F. Eisele, C. Cantrell, A. Clarke, V. Kapustin, R. 
Talbert, J. Dibb, A. Bandy and D. Thomton. A case study of the marine sulfur 
budget based on P-3 airborne observations recorded during PEM-Tropics B. 
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Appendix B 


Papers Resulting from PEM-Tropics B 

1. Nucleation in the equatorial Pacific during PEM Tropics B: Enhanced boundary 
layer H 2 SO 4 but no particle production, R.J. Weber, K. Moore, V. Kapustin, A. 
Clarke, R.L. Mauldin, E. Kosciuch, C. Cantrell, F.Eisele, B. Anderson and L. 
Thornhill. Submitted to Journal of Geophysical Research-Atmospheres. 

2. Marine latitude/altitude OH distributions: Comparison of Pacific Ocean 
observations with models, D. Davis, G. Grodzinsky, G. Chen, J. Crawford, F. 
Eisele, L. Mauldin, D. Tanner, C. Cantrell, W. Brune, D. Tan, I. Faloona, B. 
Ridley, S. Sandholm, G. Sachse, S. Vay, B. Anderson, M. Owens, B. Heikes, J. 
Snow, R. Shetter, B. Lefer, D. Blake, N. Blake, M. Carroll and Y. Wang, Journal 
of Geophysical Research, in press. 

3. Chemical ionization mass spectrometry technique for the detection of 
dimethylsufoxide and ammonia, J.B. Nowak, L.G. Huey, F.L. Eisele, D.J. Tanner, 
R.L. Mauldin III, C. Cantrell, E. Kosciuch and D.D. Davis. Submitted to Journal 
of Geophysical Research. 

4. Measurements of OH aboard the NASA P-3 during PEM-Tropics B, R.L. 

Mauldin III, F.L. Eisele, D.J. Tanner, C.A. Cantrell, E. Kosciuch, J.B. Nowak, G. 
Chen, L. Wang, D. Davis, B.A. Ridley and B. Lefer, Journal of Geophysical 
Research, in press. 

5. The relationship between OH measurements on two different NASA aircraft 
during PEM-Tropics B, F.L. Eisele, R.L. Mauldin, D.J. Tanner, C. Cantrell, E. 
Kosciuch, J.B. Nowak, B. Brune, I. Faloona, D. Tan, D.D. Davis, L. Wang and G. 
Chen. Journal of Geophysical Research, in press. 

6 . Airborne observations of DMSO, DMS, and OH at marine tropical latitudes, J .B. 
Nowak, D.D. Davis, G. Chen, F.L. Eisele, R.L. Mauldin III, D.J. Tanner, C. 
Cantrell, E. Kosciuch, A. Bandy, D. Thornton and A. Clarke, Geophysical 
Research Letters, 28, 2201-2204, 2001. 
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APPENDIX C 


Measurements of OH aboard the NASA P-3 during PEM-Tropics B 


R L Mauldin III, F.L. Eisele, D.J. Tanner, C A Cantrell, E. Kosciuch, J.B. Nowak, G. 
Chen, L. Wang, D Davis, B.A. Ridley, and B. Lefer 


Abstract. 

Airborne measurements of the hydroxyl radical were performed aboard the NASA 
P-3B using the Selected Ion Chemical Ionization technique during the Pacific Exploratory 
Mission (PEM)-Tropics B study. Typical midday boundary layer OH concentrations of 6- 
7 x 10 6 molecule cm' 3 were observed in the vicinity of Christmas Island (0-3°N), and 7-8 x 
10 6 molecule cm' 3 in the vicinity of Tahiti (17-22°S). Photochemical box model 
calculations of the concentration of OH ([OH]) throughout the entire mission yield 
generally good agreement (slope = 1.26) with a tendency for the model to slightly 
overestimate the measured OH Boundary layer measurements of OH in an air mass 
containing low ozone concentrations (<5 ppbv) yielded a midday [OH] of ~5 x 10 6 
molecule cm' 3 , some 45% lower than 7 x 10 6 molecule cm' 3 observed in air masses at the 
same latitude with more typical ozone concentrations (30-40 ppbv). Measurements 
performed under clouds yielded [OH] of ~2.5 x 10 molecule cm , a value which is a 
factor of two lower when compared to data taken on either side of the cloud. This 
decrease corresponds to the decrease in j( 03 ) to produce O(’D). 
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1. Introduction 


The hydroxyl radical (OH) is the principal oxidant in the Earth s lower 
atmosphere, reacting with a wide variety of compounds. While its concentration is 
typically less than a pptv (part per trillion by volume), its high reactivity makes it the 
primary tropospheric cleansing agent responsible for controlling the build-up of many 
compounds of both anthropogenic and biogenic origin. These compounds include such 
species as CO, CH4, numerous other non-methane hydrocarbons and their oxidation 
products. 

The primary source of tropospheric OH is the photolysis of ozone to produce 
O('D) followed by the following reaction sequence: 

O3 + hv(< 330 nm) — > O('D) + 0 2 ('A g ) j(0( 1 D)) 

(Rl) 0(‘D) + H 2 0 -4 20H 

The photolysis of H 2 0 2 and CH ? OOH can also produce OH directly but to a much smaller 
degree. Once formed, OH reacts primarily with CO and CH 4 in the remote troposphere to 
produce H0 2 and CH 3 0 2 respectively, with the latter species eventually producing H0 2 . 
The H0 2 (OH reservoir species) can react further to regenerate OH via the reactions: 

(R2) H0 2 + NO -> OH + N0 2 

(R3) H0 2 + 0 3 -> OH + 20 2 

Thus the production of OH depends primarily on the available actinic flux (which 
determine the rates of photolysis), the [0 3 ], and the [H 2 OJ; and secondarily upon the 
concentration of NO. 

Recent measurements of OH in and around clouds reveal elevated OH 
concentrations. These enhanced levels of OH are attributed to increased levels of actinic 
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flux. Model calculations show that the clear sky actinic flux could be enhanced by as 
much as a factor of two [Duynkerke and Van Weele, 1991; Madronich, 1987], While 
there has been work on the enhancement of OH concentrations above and inside clouds 
[Mauldin et al , 1997, 1998, 1999], there has been little investigation of the depression of 
OH concentrations below clouds or the impact of clouds upon column OH concentrations. 

In this work we present the results of airborne measurements of OH made aboard 
the NASA P-3B during PEM-Tropics B using the Selected Ion Chemical Ionization 
(SICIMS) technique. These results are compared to those calculated from a 
photochemical box model developed at the Georgia Institute of Technology. Other 
observations such as OH concentrations measured in the presence of low ozone 
concentrations, and above or beneath clouds are also discussed. 

2. Experimental Technique 

Measurements of OH were performed aboard the NASA P-3B using the SICIMS 
technique. A schematic diagram of the instrument used in this study can be found in 
Figure 1. The system consists of three major sections, a shrouded inlet which straightens 
and slows the air flow; an ion reaction region in which the chemical ionization reactions 
occur; and a turbo molecular pumped vacuum chamber which houses a quadrupole mass 
spectrometer and an electron multiplier detector. The measurement technique and system 
used in this study has been described in detail elsewhere [Eisele and Tanner , 1991; Eisele 
and Tanner , 1993, Tanner et al. , 1997; Mauldin et al. , 1998; Mauldin et al, 1999]; 
therefore the reader is directed to those works for experimental details. 
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The method for determining the calibration coefficient, C, in equation (I) however 
bears some discussion as a different type of calibration assembly from that used in 
previous studies by our group was employed here. Briefly, the calibration technique 
involves producing a known amount of OH in front of the 1 .9 cm curved sampling port by 
photolyzing ambient H 2 O via the 184.9 nm emission line from a filtered Hg Pen-Ray 
Lamp. In the past, a band pass filter was used to separate the 184.9 nm emission line from 
the others present. This type of filtering works quite adequately at lower altitudes where 
H 2 O concentrations are generally high. Problems arose however during the day at higher 
altitudes (>15,000 ft) where low H 2 0 concentrations allowed little OH to be produced 
from the calibration source when compared to ambient OH concentrations (except at night 
when few flights are flown). One solution to alleviate this situation would be to add a 
know quantity of H 2 0 into the sample air steam. This approach was deemed impractical 
due to the difficulties in adding large quantities of H 2 0 to large volumes of air in a known 
way. Another solution is to increase the light output of the calibration assembly. The 
bandpass filter passes 184.9 nm light with -10% efficiency. If this efficiency could be 
increased while still maintaining the same level of rejection of the other emission lines, it 
would be possible to produce a larger known amount OH under low H 2 0 conditions. 

A schematic diagram of the calibration assembly used in this study can be found in 
an insert in Figure 1 . As can be seen, light passes from a temperature controlled Hg lamp 
through a quartz cylindrical lens (focal length = 2.54 cm). The light is then reflected off 
two mirrors mounted at 45° to the light path. These mirrors have been coated to 
selectively reflect 184.9 nm. Light exits the calibration housing through a Suprasil 
window and combination shutter/slit assembly, then passes through a 2.5 cm enclosed path 
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before illuminating the sample flow. The calibration housing is N 2 purged to prevent light 
absorption, and the build up of 0 3 from the photolysis of 0 2 . Laboratory tests using this 
arrangement show an increase in the intensity of the 184.9 nm emission line in front of the 
curved inlet of a factor of 8 over the 184.9 nm bandpass filter used in previous studies. 
The shutter assembly has also been modified to allow for two photolysis settings (slit 
widths). The “narrow” setting (~1 mm) was approximately a factor of eight less intense as 
the “wide” setting (~8 mm) and mimicked the output of the calibration assembly used in 
the past. In practice it was found that at dew points greater than -10 °C, the wide 
setting produced such high concentrations of OH that the N0 3 reagent ion began to be 
depleted significantly. Thus in general, low altitude calibrations used the “narrow” setting 
while the “wide” setting was used for high altitude determinations. As discussed later, 
calibrations performed under conditions where both slit settings produced reasonable 
concentrations of OH yielded good agreement between the two settings. One additional 
feature added to the inlet system for this study is a UV sensitive photodiode mounted on 
the inner 5.1 cm flow duct opposite the calibration light entrance (Figure 1). The output 
of this photodiode is not used in the calculation of the calibration factor but rather to 
monitor the output of the calibration assembly, permitting in flight checks of the exit 
window quality and output reproducibility as a function of altitude. 

The OH concentration produced by the calibration source is a function of the 
intensity of the photon flux at 184.9 nm, the [H 2 0], the absorption cross-section at 184.9 
nm, the yield of OH from H 2 0 photolysis, and the sample flow velocity. The flow velocity 
was measured using a Pitot tube mounted just forward of the 1 .9 cm curved inlet, and the 
[H 2 0] was measured using a dew point hygrometer aboard the aircraft. The photon flux 
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at 184.9 nm was mapped out on the ground periodically during the study using a vacuum 
UV diode mounted on an x/y movable traverse. This vacuum UV diode was 
intercompared with a National Institute of Standards and Technology (NIST) standard 
diode prior to, and at the end of the mission. Recent measurements of the absorption 
cross-section of H 2 0 at 184.9 nm by Creasy et al. [2000], Hofzumahaus et al. [1997], and 
Cantrell et al. [1997] are in good agreement and yield values of 7.22x1 O' 20 , 7.0x1 0‘ 2 °, 
7. 14x1 O' 20 cm 2 molecule 1 respectively. Yoshino et al. [1996] report a value of 1 .02x10 19 
cm 2 molecule 1 , however this value was obtained near the end of their detection limit. For 
this work, the value of 7. 14x1 0‘ 2 ° cm 2 molecule 1 from Cantrell et al. [1997] was used as 
it is in good agreement with other measurements, and this value has been used for the 
reported values in earlier work [i Mauldin et al., 1998; Mauldin et al., 1999], It should be 
pointed out that all of these recent cross-section measurements are some 30-80% larger 
than the currently recommended value of 5.5x10 20 cm molecule [DeMore et al., 1997], 
a value based on the analysis by Hudson [1974] of measurements performed over 30 years 
ago. 

Calibrations were performed periodically at various altitudes during each flight. 
The photon fluence was checked after every second or third flight. With the exception of 
one measurement, measurements of the lamp output were found to be constant to within 
6% of each other for a given slit setting. The one measurement that fell outside of this 
range was not used and was caused by a dirty output window. This change was detected 
also by the monitoring diode. A calibration involved measuring the ambient OH 
concentration, then opening the shutter to a given slit setting allowing the photolysis of 
ambient H 2 0. The calculated OH produced from the calibration source was then 
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compared to the signal rise observed after the shutter was opened. Figure 2 is a plot of 
the measured calibration coefficient as a function of ambient pressure. Each point is a 
single 30 s, measurement. The scatter in the points is thought to be due to the uncertainty 
in the H 2 0 concentration. This value is used twice in the calculation of the calibration 
coefficient. First it is used to calculate of the attenuation of the lamp output reaching the 
sampling inlet. Second it is used directly in the photolysis equation used to calculate the 
[OH] produced. The dew point hygrometer aboard the aircraft had a response time on the 
order of 10’s of seconds. While not desirable, as a matter of practicality the dew point 
was found to change during some calibrations producing the observed scatter. 

The plot the points are grouped into two categories, high altitude lens setting and 
low altitude lens setting each with narrow and wide slit values. As can be seen at altitudes 
where both narrow and wide slit settings were used, the calculated calibration values are 
indistinguishable During the flights it was found that the voltages on outside atmospheric 
pressure electrostatic lenses needed to be adjusted to maintain high signal levels as the 
ambient pressure changed with changing altitude. Typically these voltages only changed 
between two different settings, one for high altitudes and one for low. While controlling 
the observed signal level, these settings also determine the lifetime in the ion reaction 
region. Thus, as a consequence of changing these settings, the efficiency that H 34 S0 4 ‘ 
(and HS0 4 ) ions are formed is also changed. This fact is revealed in the calibration values 
measured at 700 mb shown in Figure 2. The average of the high altitude lens values is 
6.00 x 10 9 while the average of the low altitude settings is 7.75 x 10 9 molecule cm 3 , 
indicating that the high altitude settings are more efficient at producing H 34 S0 4 ' ions. A 
linear fit to calibrations taken at either lens setting is also shown in Figure 2. Depending 
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upon the lens values recorded during the measurement, either the fit to the high or low 
altitude data was used to obtain the OH concentrations presented here. Unfortunately, 
these lens voltages were not recorded at all times during the mission. OH concentrations 
measured during these periods were obtained either by using the average value of the 
calibration coefficient calculated for the two different lens settings for data taken at 
altitudes with pressures less than 850 mb or by using the low altitude lens fit for data taken 
at altitudes with pressures greater than or equal to 850 mb. The total reported error limits 
of ±40.0% are 2 sigma and correspond to the 95% confidence level referring to the total 
systematic and random errors for a given measurement. The uncertainty was obtained by 
a thorough propagation of error calculation Data taken during periods where no lens 
voltages were recorded have an additional uncertainty of ±20% or a total uncertainty of 
±60%. This increased uncertainty reflects the use of a “best guess” calibration value for 
that data. 

It should be noted that high altitude OH concentrations presented in earlier work 
[Mauldin et al, 1998; Mauldin et al, 1999] were calculated using an exponential fit to 
calibration data. This fit began at a value of 9.8 at sea level, and rose to 14.8 at 465 mb. 
As stated earlier, it was difficult to perform calibrations at high altitudes during the 
daylight. As a result there was a large uncertainty in those calibration values obtained. In 
light of our present set of high altitude calibration values, the authors now feel that a 
lower calibration value should be used than was originally applied. A more appropriate 
value would be the line shown as Linear Fit A in Mauldin et al. [1998] and Mauldin et al. 
[1999], Using this fit will not change the low altitude values, as the deviation ffom the 
exponential fit does not begin until -10,000 ft. At 10,000 the correction is <4%, rising to 
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a maximum correction at altitudes greater than or equal to 18,000 of -38%. While this 
change is significant, it applies only to data taken at higher altitudes and not to boundary 
or buffer layer measurements where most of these studies were conducted. 

The results presented here are given as 30 s measurements. A typical 30 s 
measurement consisted of an 8 sec measurement of mass 99 (H 34 S0 4 ), with only 34 S0 2 
and N 2 being added through the front set of injectors in the curved sampling port. At the 
end of this 8 s measurement (called the OH signal measurement), the flow of N 2 is 
replaced by propane to remove any OH present. While waiting for this flow to equilibrate, 
mass 62 (N0 3 ), mass 97 (HS0 4 ), and an electronic background were measured for 1, 4, 
and 2 s respectively. Next mass 99 was measured again for 8 s (called the OH 
background measurement). At the end of this OH background measurement, the propane 
flow is turned off and N 2 again added. While waiting on this flow to equilibrate, mass 62 
(NOj ), mass 97 (HS0 4 ), and mass 95 (MS') were measured for 1, 4, and 2 s respectively. 
The value of HS0 4 ' in equation (I) above is the time normalized counts at mass 99 without 
propane present in the front injector minus the time normalized counts with propane 
flowing. In this manner we report [OH], [H 2 S0 4 ], and [MSA] once every 30 s, although 
it is possible to report [H 2 S0 4 ] once every 15 s. 

3. Measurements and Modeling Results 

The results of OH measurements performed during the entire PEM-Tropics B 
campaign are shown in Figure 3. The data are presented as 30 s measurements using the 
technique described above. The gaps present in the data are due to several different 
reasons. The data shown have been filtered to remove any points deemed to have 
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unacceptable auxiliary parameters, including flows, lens settings, electrical noise, or 
background signal. Data taken during times of calibrations have also been removed, as the 
concentrations of OH produced are so large that it is not possible to reliably retrieve 
ambient OH concentrations during these periods. The remaining gaps are times where OH 
was not measured to gain better precision on other species being measured by the 
instrument such as 1 Hz H 2 S0 4 data for nucleation studies, or periods where the 
instrument was taken offline due to dense precipitation or to perform diagnostics or repair 
work. For reference purposes, an altitude profile is also presented for each flight. As can 
be seen, most of the large-scale structure or variability is due to changes in altitude. The 
smaller scale features observed in many of the flights are due to changes in [H 2 0], KO3) 
and/or changes in [NO], It should be noted that the OH high frequency scatter actually 
decreases percentage wise with increasing [OH], As this is a fast OH measurement 
technique, it is possible to capture “fast” changes due to changes in [H 2 0], While there 
were fast 0 3 and NO, there was no fast H 2 0 measurement aboard the P3-B. 

The period of 04:00 to 05.00 UTC during flight 14 was flown in the dark and is a 
good example of the precision that the instrument is capable of when measuring low [OH], 
Typical 5 min averages during this time period yield values of (3-5 ± 10) x 10 
molecule cm' 3 . Under clear cloudless conditions, the average midday boundary layer 
concentration of OH was ~6-7 x 10 6 molecule cm' 3 in the vicinity of Christmas Island (0- 
3°N) and 7-8 x 10 6 molecule cm' 3 in the vicinity of Tahiti (17-22°S). 

Photochemical box model simulations of the data shown in Figure 3 were 
performed using a model developed at the Georgia Institute of Technology (Ga. Tech ). 
The model used for these simulations is described folly elsewhere here [Wang et al., this 
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issue]; thus only relevant details will be discussed. The model inputs were one-minute 
averages of the species measured aboard the aircraft that affect OH concentrations (O3, 
j(0 3 ), NO, CO, H 2 0 2 , hydrocarbons, etc). Reaction rates were based on the 1997 NASA 
JPL [DeMore et al . , 1997] and the 1997 J. Phys. Chem. Ref. Data [Atkinson et al. , 1997] 
recommendations. 

An important new aspect to the model inputs is the values for photolysis rates or j 
values. In previous works involving OH modeling, it was necessary to use either 
corrected Eppley radiometer measurements or model calculations to determine photolysis 
rates [McKeen et al., 1997; Mauldin et al., 1997, 1998; Davis et al., 1999 and references 
therein]. In this work however measured photolysis frequencies were used as model 
inputs. These values were obtained from two 2n str (combined to give 4 tc str total 
response) spectral radiometers mounted above and below the aircraft [Lefer et al., this 
issue]. These photolysis rates should allow much better estimates of OH production, 
especially for regions in and around clouds. 

A plot of measured vs modeled [OH] for flights 4-19 (the entire PEM-Tropics B 
campaign) is shown in Figure 4. Also shown is the 1:1 correlation line for reference. As 
can be seen with the exception of a few outlying points the agreement is generally good. 
A standard regression of this data yields a slope of 1 .26 with a correlation coefficient of 
0.93. This value of the slope indicates that the model has a slight general tendency to 
over-predict OH concentrations when compared to the measurements presented here. It 
should be pointed out however that this tendency is slight and that the agreement is well 
within the uncertainties of the measurements or model inputs. 
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P-3 Flight 14 was designed to investigate the HO* chemistry at various altitudes 
from midday to sunset. The flight involved flying a wall profile with four altitudes, 1000, 
4000, 10,000, and 18,000 ft. In flying a profile of this type, the each altitude can be 
sampled with the sun at different solar zenith angles. Figure 5 is a plot of the measured 
and modeled [OH] and is a good example of the typical agreement between the two. Also 
shown are j(0 3 ) and [H 2 0], two parameters that directly control the production of OH, 
and an altitude profile. As can be seen the agreement between the measurement and 
model is generally quite good, despite the [H 2 0] and j(0 3 ) changing by over a factor of 20 
and 200 respectively throughout the flight. It is interesting to note that the largest 
discrepancy of -60% occurs at 00:15 UTC at the highest altitude. However, at 01:45 
UTC when the aircraft was again at 18,000 ft, the model and measurements show good 
agreement to the 18% level. Flight 10 also demonstrates the agreement between 
measurements and model under unusual conditions. During this flight low [0 3 ] (<5 ppbv) 
were observed in the marine boundary layer. As [0 3 ] also directly controls OH 
production, measurements and model should reflect the lower concentration. Figure 6 is a 
plot of measured and modeled OH together with observed [0 3 ] and an altitude profile. As 
can be seen, the agreement is again quite good even with the [0 3 ] changing by more than a 
factor of six as the aircraft changes altitude. The observed midday [OH] of -5 x 10 6 
molecule cm" 3 is some 45% lower when compared with the value of 7x10 molecule cm 
observed during flight 8, a flight flown at the same latitude, but with [0 3 ] of -30 ppbv (a 
more typical value). 

In Figure 3 - flight 17, the data at 20:45 UTC show extremely depressed (-2.5 x 
10 6 molecule cm" 3 ) concentrations of OH when compared to values measured ten minutes 
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before or after (~6 x 10 4 * 6 molecule cm' 3 ). These lower concentrations were observed 
during a portion of the flight in which the aircraft flew under a large deck of stratoform 
clouds. Later at 22:45-23:00, elevated concentrations of OH are observed when 
compared to the times on either side of the interval. During this portion of the flight, the 
P-3 flew above a cloud layer located at -4,500 ft. Elevated concentrations of OH have 
been observed over clouds in the past by our group [Mauldin et al . , 1997, 1998, 1999]. 
In those works we concluded that the elevated OH concentrations were caused by an 
increase actinic flux. In contrast to that work, the 20.45 time period demonstrates the 
ability of a high altitude cloud to lower [OH] by attenuating the actinic flux under it. 
Figure 7 is a plot of measured and modeled OH together with j( 03 ) values determined 
from the spectral radiometer aboard the aircraft. As can be seen during this time period 
the value of j( 03 ) falls by approximately a factor of 2.5, and the [OH] drops by about a 
factor 2.4 demonstrating the correlation between OH and KO 3 ). During the 22:45-23:00 
UTC time period, the j(0 3 ) increases by -45%, and the OH increases from -5 x 10 6 to 
-9.5 x 10 6 molecule cm' 3 . Presumably these increased values arise from increased 
upwelling radiation from the cloud top combined with slightly increased H 2 0 above the 
cloud. The modeled OH again shows good agreement, even during these periods of 
increased and decreased actinic flux. 

4. Conclusions 

The SICIMS technique described here allowed a large number of OH, H 2 S0 4 , and 

MSA measurements to be performed aboard the NASA P-3B during the PEM-Tropics B 
study. These measurements are further enhanced by the use of an improved calibration 
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assembly, which allows reliable calibrations to be performed under low [H 2 0] (high 
altitudes). The use of this improved calibration system has reduced the reported 
uncertainty of high altitude OH values by a total of 20% (42% vs 62%) compared to 
values previously reported by our group [Mauldin et al, 1997, 1998], Additionally, in 
light of the high altitude calibration values obtained in this study, the high altitude values 
previously reported by Mauldin et al. [1997, 1998] should be lowered by as much as 38/o 
for data obtained at the highest altitudes. These corrections do not apply to data 
presented here. If corrected values are desired, the reader is directed to contact the 
authors. 

When combined with other parameters measured aboard the aircraft, the OH 
measurements presented here provide an extensive data set for comparison with model 
calculations. Overall comparison with the Ga. Tech, photochemical box model 
simulations yields generally good agreement. The calculated slope of 1 .2 shown in Figure 
4 indicates that the model has a tendency to predict slightly higher values than measured. 
As shown in Figure 5, the largest discrepancies of -60% were observed at the highest 
altitudes where the [H 2 0] generally dropped to less than 1 mb. While this concentration 
of H 2 0 is not extremely small, it is nearing the detection limit for the dew point 
hygrometer aboard the aircraft. If there are errors in the reported [H 2 0], a key compound 
in the production of OH, there will be errors in the [OH] calculated by the model. The 
above argument is not meant to explain these discrepancies, but to point out the need for 
accurate measurements of species important in the production of OH when trying to 
model OH or HO x photochemistry. 
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Another highlight of the study was an encounter of an air mass with low O3 
concentrations. During flight 10, boundary layer [0 3 ] of <5 ppbv were observed as 
compared to more typical values of 30-40 ppbv. Such a suppression of a compound 
important in the production of OH is revealed in both measurements and model 
simulations. The data shown in Figure 6 reveal a midday boundary layer [OH] of ~5 x 10 
molecule cm' 3 , a value some 45% lower than the value of 7 x 10 6 molecule cm' 3 observed 
during flights flown at the same latitude, but with more typical 0 3 concentrations (flight 
8). It is interesting to note that while [0 3 ] is a factor of 6-8 lower than typical levels, the 
midday [OH] only dropped by a factor of 2. As the production of OH has a square root 
dependence on ozone concentration, the production of OH is lowered by a factor of ~2-3, 
which is in the range of the reduction of OH. Model simulations reveal that the OH 
production terms for boundary layer portions of flight 10 are approximately a factor of 2- 
2.5 smaller than calculated for flight 8. Thus, while the ozone concentration of the air 
mass encountered during flight 10 has been reduced significantly, the production level of 
OH is still maintained at levels to sustain a reasonable oxidizing capacity. 

The portion of flight 17 flown around clouds is also of interest. In addition to 
elevated levels of OH observed in and over clouds previously reported [Mauldin et al, 
1997, 1998], here we also report depressed values of OH under a cloud. In Figure 7, it 
can be seen that [OH] falls from ~6 x 10 6 molecule cm' 3 observed on either side of the 
cloud encountered at 20:45 to 3 x 10 6 molecule cm' 3 observed under it. This decrease in 
[OH] is due to a decrease in actinic flux, as can be seen by the large change in j(0 3 ) shown 
in Figure 7. For the cloud encountered at 22:45-23:00 UTC time period, the OH 
increased by nearly a factor of two as the aircraft flew above the cloud. Again this change 
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can be attributed to an increase in actinic flux as indicated by the increase in XO 3 ) 
observed during this time period. Thus, in this single flight the ability of clouds to both 
increase and decrease OH concentrations is demonstrated. While these changes of OH 
around clouds are not unexpected, they do raise the question of how clouds impact the 
column concentration of OH. In light of the present data, it would appear that cloud 
altitude plays a significant roll in such a calculation. If a cloud were at a lower altitude 
such as the boundary layer where H 2 O concentrations are typically high, it is possible that 
the reduced concentration of OH produced beneath the cloud is compensated for by 
enhanced OH production above and inside the cloud leading to either a constant or larger 
column concentration If however, the cloud were at a high altitude where H 2 0 
concentrations are typically low, it is possible that the increased OH production above and 
inside the cloud could be overwhelmed by the reduction in OH production below the 
cloud leading to a decrease in column [OH], Thus, depending upon a clouds altitude, it 
could lead to an increase, decrease, or a redistribution of column OH. Further discussion 
of the impact of clouds upon OH concentrations will be given in a future publication. 

It is also interesting to note that the modeled OH concentrations show generally 
good agreement, even during the periods of increased and decreased actinic flux. This 
agreement is much better than the simulations presented in Mauldin et al. [1997, 1998] 
for conditions in and above clouds using either the Ga. Tech model or one developed at 
the NOAA Aeronomy Lab. These models used as inputs actinic fluxes derived from 
Eppley radiometer measurements. As stated in Mauldin et al. [1987, 1998] actinic fluxes 
derived in this manner tend to under estimate j(C> 3 ) values. In contrast, the simulations 
presented here used j values obtained from the spectral radiometer flown aboard the 
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NASA P-3B. In comparison to the Eppely, which can only measure intensity in a limited 
wavelength range, this instrument can obtain direct determinations of j values at different 
photolysis frequencies, an ability that greatly enhances the model derived OH values. 
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Figure Captions 


Figure 1. Schematic diagram of the SICIMS instrument used for the airborne 
measurement of OH, H 2 S0 4 , and MSA during the PEM-Tropics B study. The inset 
shows the new calibration source, which allows better calibrations to be performed 
under low H 2 0 concentrations (high altitudes). 

Figure 2. Plot of calibration coefficient versus pressure altitude obtained during the 
PEM-Tropics B study. Measurements and calibrations were performed under two 
different electro-static lens settings. As can be seen with the 700 mb values, these 
settings affect the obtained value for the calibration coefficient. Also shown are linear 
fits to the values for each lens setting. Depending on the lens setting recorded, the 
appropriate fit was used to obtain the OH concentrations presented here (see text). 


Figure 3. Concentrations of OH measured during each flight during the PEM- 
Tropics b study. Also shown for reference is an altitude profile for each flight. Data 
are presented as 30 s measurements, which are the result of 8 s of OH signal 
measurement, and 8 s of OH background measurement during each 30 s period (see 
text). 

Figure 4. Comparison of measured OH values from the entire PEM-Tropics B study 
with those obtained using the Ga. Tech, photochemical box model. Also shown is the 
line indicating 1:1 correlation between the values. A linear fit yields a slope of 1.26, 
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indicating a slight tendency for the model to over-predict OH concentrations when 
compared to the measured values presented here. 

Figure 5. Plot of measured and modeled [OH] for P-3 flight 14. Also shown for 
reference are an altitude profile, and the values of KO3) and [H2O]. The agreement 
between modeled values and measurements is generally good with the largest 
deviations occurring at high altitudes. The last portion of the flight was flown at night. 
Note the precision in measuring low OH concentrations. Five-minute averages during 
this time period yield values of (3-5 ± 10) x 10 4 molecule cm 3 . 

Figure 6. Plot of measured and modeled OH for P-3 flight 10. Also shown for 
reference are an altitude profile and the measured [O 3 ]. Note the low ozone 
concentrations observed in the boundary/buffer layer runs. Model simulations show 
that the OH production in these regions are ~2-2.5 times smaller than those calculated 
for flight 8, a flight flown at approximately the same latitude, but in which more typical 
[O3] were observed. 

Figure 7. Plot of measured and modeled OH for P-3 flight 17. Also shown for 
reference are an altitude profile and the measured j( 03 ). Time periods where the 
aircraft encountered clouds are marked. Note the large decrease in OH observed as the 
aircraft passed beneath the cloud at 20:45. A large increase in OH is observed as the 
aircraft flew above a cloud at 22:45-23:00. This flight demonstrates the ability of 
clouds to both increase and decrease OH concentrations. Also note the good 
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agreement of measured and modeled OH values, even during the periods of increased 
and decreased actinic flux. 
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APPENDIX D 


The Relationship Between OH Measurements on Two Different NASA Aircraft 

During PEM-Tropics B 


F. L. Eisele, R. L. Mauldin, D. J. Tanner, C. Cantrell, E. Kosciuch, J. B. Nowak, 
B. Brune, I. Faloona, D. Tan, D. D. Davis, L. Wang, and G. Chen 


Abstract. OH measurements were performed on both the P3-B and the DC-8 aircraft 
throughout PEM-Tropics B, using two very different measurement techniques. While no 
direct comparison was possible because of the difference in flight paths of the two 
aircraft, 2 brief measurement periods in close proximity allow at least a glimpse of how 
the two measurements compare. The first comparison took place in the marine boundary 
layer and showed exceptionally good agreement between the two aircraft OH 
measurements. The second set of close proximity flights at 5.5 km altitude resulted in an 
average concentration difference approximately equal to the uncertainty limit associated 
with either individual instrument, but much less than the combined uncertainties from all 
the contributing measurements. In both cases the comparisons were made by 
normalizing the data with a photostationary state model. In addition, a comparison of the 
OH concentrations measured by each aircraft from 25 degrees north latitude to 25 
degrees south latitude over the entire longitude and altitude range of the mission resulted 
in agreement within about 10% but with a trend of higher DC-8 /P-3B OH ratios at higher 
altitudes. While a definitive comparison of these two OH instruments must await a far 
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more rigorous future test, the present results show that there were no clear OH 
measurement discrepancies observed. The dramatic differences in these two 
measurement techniques also minimizes the chance of common interferences or 
calibration errors. 


1. Introduction 


The hydroxyl radical is the primary cleansing agent of the lower atmosphere. It is 
initially derived from the photolysis of ozone to form 0 ( l D) and its subsequent reaction 
with water to yield two OH radicals. Once formed, OH reacts with most compounds 
emitted into the atmosphere, including virtually all hydrocarbons and their oxidation 
products and most sulfur and nitrogen containing compounds. Depending on the specific 
reaction, OH will either be destroyed as, for example, when it reacts with NO2 to form 
HNO3, or recycled, as when it reacts with CO or methane to form HO2 or ROi. The 
highly reactive nature of OH makes it difficult to measure and to model for several 
reasons. In the case of measurements it can react with relatively stable and abundant 
gases such as CO and methane, thus its lifetime is short, on the order of a second, and its 
peak concentrations are typically held very low, i.e., in the 10<> molecules cm*' range. 
Measurements must therefore be accomplished quickly, sensitively, and without 
significant new reaction possibilities such as contact with inlet or sample tube walls. On 
the modeling side, the primary source of OH is fairly well established, at least at low 
altitudes, but the sinks include CO, methane, plus small amounts of numerous other 
compounds. Many of the reaction rates with these other compounds are unknown, their 
concentrations are often not measured, and their reaction products are also often in 
question. There are also several secondary sources of OH which involve various forms 
of chemical and photolytic recycling through compounds such as HO2, H2O2 and their 
organic counterparts. 
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Despite measurement difficulties, it is critical that OH be accurately measured 
because it is the concentration of this radical that directly controls the lifetime of most 
biogenic and anthropogenic emissions. It is also important that models accurately portray 
the role of various parameters in controlling OH concentrations if the consequences of 
changes in anthrogenic or biogenic emissions are to be properly understood. Before 
model/measurement comparisons can be meaningfully discussed, however, it must first 
be determined how models compare with each other (which is beyond the scope of this 
paper) and how measurements on the two aircraft compare. It is the purpose of this paper 
to provide some new insight into the level of agreement between the two aircraft OH 
measurement instruments used in PEM-Tropics B. The authors will be the first to admit 
that the data to be discussed does not constitute a rigorous intercomparison. It is, 
however, important to use the limited data set in which OH was measured on both aircraft 
at similar locations and times to make sure that agreement between the two very different 
OH measurement techniques is at least sufficient to allow meaningful comparisons of the 
photochemical results from the NASA DC-8 and P-3B aircraft. If the two OH 
instruments consistently differ by more than their combined limits of uncertainty when 
measuring under apparently similar conditions, then clearly the interpretation of 
measurement/model comparisons would be very different, depending on which aircraft 
platform the data came from. If, on the other hand, measurements typically agree within 
the limits of uncertainty of each instrument, then there is at least no a priori expectation 
that measurement/model interpretations are platform dependent. This does not mean that 
the two instruments can be said to agree or even that they are not both wrong. A 
statement of such agreement or lack thereof awaits a far more rigorous intercomparison. 
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The possibility of both instruments being biased in the same way or containing similar 
interferences may be even more difficult to completely rule out. These two techniques, 
however, are so different in operating pressure, analytical methods, and ion or photon 
detection that it would be difficult to conceive of a mechanism or chemical interference 
that might be responsible for a common bias in both instruments. 

2. Instrumentation and Measurement Techniques 

The instrumentation used to measure OH on the DC-8 and P-3B could not be 
much more dissimilar. The DC-8 instrument uses laser induced fluorescence (LIF) to 
preferentially excite OH molecules and then employs photon counting to detect 
fluorescent photons emitted by the excited OH molecules after a brief time delay. The 
measurements take place at a reduced pressure of only a few torr, which is achieved by 
rapid pumping of sampled air through a critical orifice. The aircraft inlet design is unique 
to this instrument and is well-modeled. Calibration is achieved by generating a known 
concentration of OH via several laboratory techniques including H 2 O photolysis at 184.9 
nm and measuring them at a range of internal pressures from about 0.3 to 1.0 kp. This 
pressure range corresponds to that in the measurement cell at maximum flight altitude 
and sea level respectively.. Relative changes in detection sensitivity are also monitored 
using 184.9 nm photolysis of H 2 O in flight. The DC-8 instrument is located in the 
forward belly of the aircraft and samples air directly below the fuselage. Measurements 
are typically made at mach 0.7± 0.1 and altitudes of 15,000 feet to 41,000 feet, with 
profiles down to the ground 
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The P-3B instrument uses chemical titration of ambient OH with ^S02 to produce 
H2^ 4 S04 which is subsequently measured by selected ion chemical ionization mass 

spectrometry (SICIMS) using NO3' reactant ions. The titration and chemical ionization 
takes place essentially at ambient pressure and temperature, and only after the ion 
chemistry is completed are the ions sampled through a N2 sheath gas into a high vacuum 
environment where they are mass analyzed. Detection is by single ion counting and OH 

concentrations are derived from the ratio of H^SO^ and NO3 ion count rates which 
leads to a very stable calibration coefficient. Calibration is achieved in situ in ambient air 
by photolyzing a known concentration of H2O in the inlet at a known velocity with a 
known photon flux at 184.9 nm photons. The inlet design used by this instrument is also 
unique to this application and has undergone extensive wind tunnel testing. The P-3B 
instrument samples near the front of the aircraft on the port side at window height. 
Measurements are typically made at mach 0.3+0. 1 at altitudes of jOO feet to 25,000 feet 
with a focus on boundary layer studies. Far more detailed descriptions of the DC-8 LIF 
instrument [Stevens et al, 1994; Mather et al., \991,Brune etal., 1998, Tanetal., 1998, 
and Faloona et al., 2000] and the P-3B SICIMS instrument [Eisele and Tanner, 1991, 
Tanner and Eisele, 1995; Eisele et al., 1996; Eisele et al., 1997; Mount et al., 1997] have 
been published previously, and the PEM-Tropics B results from this technique are 
included in this special issue. 

The similarities between these two instruments are few. They both have fairly 
similar detection sensitivities and in each case this is achieved by counting discrete 
events: i.e., the detection of individual photons or ions. Absolute measurement 
uncertainties are similar (about ± 40%), as is the use of 184.9 nm photolysis of H 2 0 as a 
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calibration standard. Both instruments have also participated in several previous aircraft 
and ground based studies [Eisele and Tanner, 1991 ,Eisele etal . , 1994, 1996; Tanner et 
al., 1997; Jefferson etal., 1998; Mauldin et al., 1998, 1999; Mather et al., \991,Bruneet 
ai, 1998, Tan et al., 1998, and Faloona et al., 2000], Apart from these functional 
similarities, the two instruments are profoundly different. The DC-8 instrument uses 
electronic excitation of OH via photo absorption at 308 nm and the detection of the 
subsequent resonance fluorescence to detect OH molecules. The P3-B alternatively 
makes use of the chemical reactivity of OH with SO 2 to form a product which is 
preferentially detected because of its high acidity and specific molecular mass. Thus, in 
order for these two instruments to have a common interfering species that could be 
mistaken for OH in each, this compound would need to have an electronic excitation 
level similar to OH, form a high acidic compound with 34 S0 2 , and result in an ton with 

the same mass as H 34 SO;. Other differences include operating pressures. The DC-8 
instrument samples air through an orifice into a low pressure region (about 0.4-1. 1 kp) 
where the LIF detection if accomplished. On the P-3B, the chemical reactions used to 
detect OH are all accomplished at pressures slightly above ambient due to ram air 
pressure in the inlet. The resulting ions are sampled into the low pressure mass analysis 
region only after the reactions are stopped by passing the ions through a N 2 buffer gas 
just in front of the vacuum entrance aperture. The transition from ambient to low 
pressure is therefore also quite different, in one case being in air and involving the OH 

radical directly and in the second being in dry N 2 and involving NOj and HSO} 10 ns. 
Again, the opportunity for water vapor, temperature, pressure or some trace gas to affect 
both of these compounds in the same way during the pressure reduction phase is not 


likely. The two instruments each have unique inlets, however, they both slow the 
sampled air flow significantly compared to free stream flow. The DC-8 operates at 
speeds up to about twice those of the P-3B, and as a result DC-8 ram air pressure and air 
temperature rises are about a factor of 3 to 4 times those in the P-3B inlet. As mentioned 
in the introduction, calibration methodology is also quite different, though both 
instruments do use 184.9 nm photolysis of H 2 O for determining a calibration factor. 

3. Measurement Results 

PEM-Tropics B was not intended to provide an intercomparison of measurement 
techniques on the DC-8 and P-3B aircraft. The two aircraft were intended to complement 
each other and not to fly redundant measurement patterns. As a result, there were few 
instances, where the two aircraft flew in a similar location at a similar time and probably 
none in which they actually flew through the same air mass. There were, however, two 
instances on flight #6 on the P-3B and flight #8 on the DC-8 and on flight #14 on the 
P-3B and flight #16 on the DC-8) in which similar air masses could have been sampled. 
The first potential overlap was near Christmas Island, when the DC-8 flew at the same 
altitude through a similar geographical location to that sampled an hour earlier by the P- 
3B. Table 1 shows a comparison of measurement times, locations, photolysis rates and 
the average concentrations of several compounds associated with either OH production or 
loss measured on the P-3B and DC-8. Most of the compounds measured on the two 
aircraft look quite similar except NO, which is on average 4 times higher on the P-3B 
than the DC-8. NO measurements on both aircraft, however, involved concentration 
spikes that were also reflected in the OH concentration. This suggests that while the air 
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masses are fairly similar in several respects they are not really the same. More than a 
factor of 2 difference in average JO('D) (but much less in J(N 02 )) also suggests that a 
direct comparison of measurements from the two aircraft is not possible. Rather than 
trying to correct for differences in individual parameters, a model will be used as a 
transfer standard. The model is a photochemical box model containing full HOx-NOx- 
CH 4 chemistry and parameterized NMHC chemistry. This model has been previously 
described in detail by Crawford et al. [1999], Instead of comparing measured OH 
concentrations directly, the ratio of measured to modeled OH on each aircraft is 
compared. Since the measurement conditions were quite similar for the two aircraft, the 
use of a model as a transfer standard should introduce little additional uncertainty due to 
uncertainties in the reaction and photolysis rates used. The uncertainties associated with 
the measurement of the most important photochemical parameters on each aircraft do, 
however, directly contribute to the uncertainty of the comparison. Table 1 also shows the 
average measured OH concentrations and the ratio of measured to modeled OH for each 
aircraft. Figure 1 shows the temporal dependence of the data sets from both aircraft. 

Note that there are spikes in both data sets, making comparison even less certain. The 
very high level of agreement of average measurement/model OH ratios shown in Table 1 
has little significance since calibration uncertainties are much larger. Clearly, however, 
agreement is very good and it is comforting to see that both instruments and models 
appear to consistently track each other well within the combined uncertainty limits of the 
measurements and model inputs. In this boundary layer comparison, both measured OH 
concentrations seem to indicate that modeled values are consistently low by about 30%. 
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This has only limited significance, since it is within the uncertainty of the OH 
measurements. 

A second opportunity for comparing DC -8 and P-3B OH measurements occurred 
on a flight out of Tahiti at an altitude of about 5.5 km. In this case, the P-3B flew 
through an exceptionally uniform air mass 5 times from about noon until late afternoon, 
local time. The DC -8 flew through nearly the same location and at the same altitude in 
the mid-afternoon. A comparison of time, location and measured parameters is given in 
Table 2 . While the mixing ratios of O 3 , CO and H 2 O were fairly consistent between both 
successive P-3B measurements and the DC -8 measurements, average NO was again 
higher on the P-3B, this time by a factor of about 2. CH 2 O and H 2 O 2 also varied by up to 
a factor of 2 on both P-3B sampling legs. This suggests that the two aircraft were 
sampling somewhat similar but still not the same air mass. Figure 2 also shows the time 
sequence of OH measurements made by the two aircraft. Several of the P-3B and the 
DC -8 data show significant structure in both the measurement and modeled OH 
concentrations, again making comparisons more difficult. At this 5.5 km altitude 
comparison, measured-to-modeled OH ratios were quite different for the P-jB and DC- 8 . 
An average of the first three or four P-3B measurement/model ratios is about 0.65, 
compared to a ratio of 1 .05 on the DC- 8 . Thus, the P-3B ratio was 62% of the DC -8 
ratio. The uncertainty in the OH measurements is about 40% for either OH instrument 
and the other measurements used in model calculation on the two different aircraft 
probably add an additional 1 5-20% uncertainty per platform to the comparison 
uncertainty. Much of the latter uncertainty comes from uncertainties in NO 
concentrations which were measured close to the detection limit. Compounds with 
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longer lifetimes such as CO and O 3 can generally be measured more accurately and 
appear to be quite similar when measured by both aircraft. The actinic flux 
measurements also have relatively small quoted uncertainties and measurements were 
made on both aircraft with essentially identical instruments. Therefore the observed 
disagreement in OH between the two aircraft is well within the combined instrumental 
uncertainty. It is even within or very close to the uncertainty imposed by the sum of the 
errors from either platform individually. Thus, the present result does not suggest a 
discrepancy between these two instruments. Neither, however, does the data suggest 
consistent agreement. With variations in both measured and modeled OH of about a 
factor of 2 within several individual flight tracks, and with variations of a factor of 2 in 
the average concentration of short-lived compounds such as NO between the two aircraft, 
it is not at all clear that the same air mass is being sampled. Thus, the comparisons of the 
ratio of measured/modeled OH for similar conditions (as determined by other measured 
photochemical parameters) provides additional evidence that the two instruments are 
measuring the same compound and that at least in these two situations they agree well 
within their combined uncertainty limits. These results are also consistent with previous 
ground-based intercomparison results for these two techniques [Mount et al., 1997 and 
Hofzumahaus et al, 1998], For these two aircraft comparison periods, however, it is not 
clear that only close spatial and temporal proximity of the aircraft added significant new 
insight into how well the measurement techniques compare. Unless the two aircraft fly 
side by side for 10 or 20 minutes, and are truly co-located it is unlikely that proximity 
provides significant benefits. Since this was not the case in PEM-Tropics B, another 
alternative is to try to compare the bulk measured OH data ratio over a very large area. 
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Davis et al., this issue, have sorted the OH data from the P-3B and DC-8 into altitude and 
latitude bins. They have used relatively clear sky measurements and adjusted them for 
diurnal changes in actinic flux to produce a map of OH concentrations expected at local 
noon in each of these bins. Averaging the median data from each latitude block common 
to both aircraft missions (which includes all longitudes) and dividing DC-8 diumally 
adjusted OH data by P-3B diumally adjusted OH data produces the altitude dependent 
OH measurement ratio shown in Figure 3. While the OH ratios largely stay within the 
uncertainty bounds of either instrument, a clear altitude trend is observed. The 
significance of Figure 3, however, should be evaluated in light of several issues. The 
highest OH ratio corresponds to a very limited latitude range at the highest P-3B altitude. 
The DC-8 measurements extended much further to both the east and west than did those 
of the P-3B. Finally, the altitude dependence of OH measurements on each aircraft may 
contain different biases. The DC-8 flies at low altitudes typically as part of a vertical 
profile in a selected area. The P-3B spends most of its time involved in low altitude 
studies and typically chooses high altitude flight paths to optimize transit time. 

4. Conclusions 

Intercomparing instruments on two very different aircraft can pose substantial 
difficulties. Differences in air speeds, operating altitudes, modification of airflow by the 
air frame, and the dangers of flying too close together or in the wake of another aircraft 
all add to these problems. The large number of diverse instrument types operated on the 
two aircraft by many different investigators, makes comparative data interpretation even 
more susceptible to error. This paper provides some insight into the relation between OH 
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measurements on the P-3B and DC-8 aircraft but a far better comparison is needed. 

Fairly consistent agreement well with the combined uncertainties of the two instruments 
provide some evidence that measurement/model comparison on both aircraft are 
providing useful new insights. DC-8/P-3B OH ratios of about 1.5 at higher altitudes both 
in adjacent flight paths and in the averages shown in Figure 3, however, suggests the 
need for a better means of intercomparison at several different altitudes in the future. 

After looking at the two sets of adjacent flight track data it also appeared that a 
comparison of all short and intermediate lifetime species could be highly enlightening. It 
is clear, however, that unless the two aircraft fly much closer together in space and time 
(maybe side-by-side within 1-2 km at altitudes where this is possible), that general 
proximity offers little assurance of flying in the same air mass. To critically compare 
multi-aircraft measurements there needs to be some assurance that the same air mass is 
being sampled and this requires more than just similar concentrations of long lived 
compounds and j values. Flights sufficiently close that they both observe spikes or 
perturbations with similar relative concentration changes, and spatial extent would, 
however, help provide evidence for a similar sampling environment. Direct comparison 
of similar measurements on the same aircraft would, of course, provide the best means of 
comparing instrument operation. 
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Figure Captions 


Figure 1 . OH concentrations measured by the P-3B on flight 6 and the DC-8 on flight 8 
as a function of time. The large differences in concentration are mainly due to 
differences in J values and NO. The OH concentrations shown are 60 second averages. 

Figure 2. OH concentrations measured by the P-3B on three legs of flight 14 and one leg 
of the DC-8 flight 16 plotted as a function of time. The OH concentrations shown are 60 
second averages. 

Figure 3. The ratio of noontime equivalent OH measured by the DC-8 divided by that 
measured on the P-3B over the P-3B altitude range. 

Table Captions 

Table 1. Comparison of average measurements on P3-B flight and DC-8 flight 8 during 
flight legs in a similar location. 

Table 2. Comparison of average measurements on several legs P3-B flight 14 and the 
only leg DC-8 flight 16 in a similar location. 
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APPENDIX E 


Chemical Ionization Mass Spectrometry Technique for the Detection of Dimethylsufoxide 
and Ammonia 

J.B. Nowak 1 , L.G. Huey 1 , F.L. Eisele u , D.J. Tanner 21 , R.L. Mauldin III 2 , C. Cantrell 2 , E. 
Kosciuch 2 , and D.D. Davis 1 

'School of Earth and Atmospheric Sciences, Georgia Institute of Technology, Atlanta, GA, 
30332, USA 

Atmospheric Chemistry Division, National Center of Atmospheric Research, Boulder, CO, 
80303, USA 

Abstract- A chemical ionization mass spectrometer (CIMS) was used to study reactions of 
protonated ethanol clusters (C 2 H 50 H) c H with dimethylsulfoxide (DMSO), dimethyl sulfone 
(DMS0 2 ), ammonia (NH 3 ), and a series of nonmethane hydrocarbons (NMHCs) and volatile 
organic compounds (VOCs). The reactivity of the (C 2 H 5 OH) n H~ cluster ions is found to be a 
function of cluster size with reactivity decreasing as cluster size increases. Ethanol cluster ion 
distributions formed from 24 ppbv, 900 ppmv, and 1 % ethanol/N 2 gas mixtures were studied. 
Small (C,H 5 OH)„H' clusters, those formed using the 24 ppbv ethanol/N, mixture, react at or near 
the collisional rate with DMSO, NH 3 , acetone, and MVK. The effective ion molecule rate 
coefficients are: 1.8 x 10* 1.5 x 10*, 1.0 x 10* and 1.6 x 10* cm 3 molec' 1 s’ 1 , respectively. 

Only DMSO and NH 3 react with the two larger (C 2 H 5 OH)„H 4 cluster ion distributions studied. 
The effective rate coefficients for DMSO and NH 3 with the 900 ppmv ethanol cluster ion 
distribution are 1.5 x 10* and 0.7 x 10* cm 3 molec’ 1 s’ 1 , respectively. The effective rate 
coefficient for DMSO with the 1% ethanol/N 2 mixture is 0.35 x 10* cm 3 molec 1 s*\ while NH 3 
reaches an equilibrium with this cluster ion distribution. Experiments show that large 
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(CjHjOH^H* ion clusters must be used at relative humidities greater than 50 % at 20 °C to 
prevent formation of and subsequent interferences from H 3 0 + ions. These results demonstrate 
that the (CjHjOH^FT ion chemistry can selectively detect DMSO and NH 3 under most ambient 
conditions with high sensitivity. 
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1. 0 Introduction 


Chemical Ionization Mass Spectrometry (CIMS) has been used to selectively measure 
many atmospheric trace gases with high sensitivity. The list of species includes: the hydroxyl 
radical (OH), sulfuric acid (H 2 S0 4 ), methane sulfonic acid (MSA), and nitric acid (HN0 3 ) 

[Eisele and Tanner, 1991 \Hueyetal, 1998; Mauldin etal., 1998; Viggiano, 1993], One 
important feature of the CIMS method is its ability to perform measurements on fast time scales 
(Le., 1 sec). For many trace gases fast time resolution measurements are essential for testing our 
understanding of sources and sinks. One such species falling into the above category is 
dimethylsulfoxide (DMSO). 

DMSO is believed to be one of the major products of the OH-initiated oxidation of DMS 
[Arseneetal., 1999; Davis et al, 1999; Davis et al, 1998; Tumipseed et al, 1996], 

OH + DMS S0 2 MSA, DMSO, DMS0 2 , other products (1) 

DMS is a by-product of marine phytoplankton metabolic activity and its release from the ocean 
surface accounts for the largest fraction of biogenic gas phase sulfur emissions [Bates et al., 

1992; Berresheim et al. , 1995]. It has been hypothesized that the formation of sulfate aerosols 
from DMS oxidation products may significantly influence the Earth’s radiation budget and 
possibly regulate climate [Charlson et al., 1987]. DMS oxidation occurs by a complex 
mechanism with multiple steps that are currently the focus of many laboratory and field 
investigations [Davis et al., 1999; Yin et al., 1990]. Critical to assessing the affect of DMS 
oxidation on climate is understanding the chemical and physical processes that control the fete of 
these products in the atmosphere. 

In the atmosphere DMSO reacts rapidly with OH 
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DMSO + OH products (2) 


There is, however, considerable debate as to the products of reaction 2 under atmospheric 
conditions. Barnes et al. [1989] and Sorensen el al. [1996] suggest that reaction 2 produces S0 2 
and DMSOj, though they disagree as to the yield of each product, while Urbanski et al. [1998] 
suggest that MSIA is a significant product. Though there is disagreement on the product yields 
of reaction 2, all the studies agree that reaction 2 leads to a DMSO photochemical lifetime of less 
than one hour for typical midday marine boundary layer OH levels. Another loss process for 
atmospheric DMSO is scavenging by cloud droplets or sea salt aerosol [Davis et al., 1999; Davis 
et al., 1998; Debruyn et al, 1995]. Therefore, depending on aerosol loading, the overall lifetime 
of DMSO in the marine boundary layer can be even shorter than that due to photochemistry 
alone. The short photochemical lifetime highlights the need for fast time resolution 
measurements of DMSO in order to gain a comprehensive understanding of DMS oxidation. 

There are few measurements of gas phase DMSO reported in the literature. The 
techniques utilized include: gas chromatography, gas chromatography/mass spectrometry, and a 
Cofer mist chamber [Bandy et al, 1996; Harvey and Lang, 1986; Sciare and Mihaiopoulos, 
2000]. One common characteristic of all of these methods is a sampling resolution of 30 minutes 
or longer. Under these conditions a detailed study of DMS oxidation becomes quite difficult. 

[Berresheim et al., 1993] developed a CIMS method to measure DMSO on a sub-minute 
time scale using NH/ core reagent ions. NH 3 has a proton affinity (PA) (204 kcal/mol) higher 
than water ( 1 66 kcal/mol) and lower than DMSO (2 1 1 kcal/mol) [Hunter and Lias, 2000], 
Therefore, NH 4 + ions proton transfer with DMSO in ambient air to form DMSOH + . The 
estimated limit of detection for a 60 s integration period has been reported at 0.5 pptv. A 
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limitation of this technique, however, is the inability to measure NH 3 . 

NHj is an another important atmospheric trace gas for which there are few fast time 
response methods available. The measurement techniques include: photofragmentation/laser- 
induced fluorescence, denuder sampling/chemiluminescence detection, denuder sampling/ion 
chromatography, and filter pack sampling/colorimetric analysis [Williams ei al., 1992]. Recent 
laboratory and field measurements have led to the hypothesis that ammonia is potentially 
important in moderating particle nucleation and growth [Coffman and Hegg, 1995; Weber et al . , 

1 998]. The primary goal of this work has been to find an appropriate reagent ion for detecting 
both DMSO and NH 3 using the same technique. One possibility explored here is protonated 
ethanol, (C 2 H 5 OH)H\ 

Ethanol has a proton affinity ( 1 85.6 kcal/mol) higher than water and lower than both 
DMSO and NH 3 [Hunter and Lias, 2000], However, (C 2 H 5 OH)H + efficiently clusters with any 
available ethanol to form an equilibrium distribution of protonated ethanol cluster ions, 
(C 2 H 5 OH) n H\ where the range of n depends on temperature and the concentration of C 2 H 5 OH. 
However, the reactivity of (C 2 H 5 OH) n H + is most likely quite from that of (C 2 H 5 OH)H\ e.g., 
dependant on, n, the cluster size [Feng and Lifshitz, 1995]. 

As a starting point in developing a CIMS method for detecting both DMSO and NH 3 , the 
effective rate coefficients for reactions of (C 2 H 5 OH) n H + with DMSO, NH 3 , DMS0 2 , DMS, a 
series of non-methane hydrocarbons (NMHCs), and a series of volatile organic compounds 
(VOCs) were measured (see Tables 1 and 2). The sensitivity of this ion chemistry for the 
detection of DMSO and NH 3 and the possibility of artifact DMSO production from ambient 
DMS oxidation within the instrument were also examined. Finally, a new ion-molecule reactor. 
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specifically designed for aircraft applications, was characterized and then used to measure 
DMSO in ambient air on the NASA P-3B as part of the PEM-Tropics B field program. 

2.0 Experimental 

2.1 Gas Standards 

Isotopically labeled C 13 DMSO and C 13 DMS0 2 permeation devices were manufactured 
by VICI Metronics, Inc. using 99 % pure isotopically labeled compounds from Icon Services. 
These devices were maintained at a constant temperature and pressure. The concentrations of 
DMSO and DMS0 2 were controlled with a dynamic dilution system similar to that described by 
[Goldan et al., 1986], The permeation rates of each device were determined gravimetrically by 
the manufacturer and are traceable to NIST standards. The combined uncertainty of the C 
DMSO permeation tube and dilution system is estimated at ± 50% and primarily reflects the 
uncertainty associated with the emission rate of the permeation tubes. The combined uncertainty 
of the C 13 DMS0 2 permeation tube and dilution is estimated at ± a factor of 5 and reflects 
additional difficulties associated with that particular permeation device. 

NH 3 concentrations were determined by monitoring the absorption of the 185 nm 
mercury atomic line. An NH 3 absorption cross section at 185 nm of 4.3 x 10 cm molecule 
was used [Lovejoy, 1999]. The nonmethane hydrocarbons (NMHCs) and volatile organic 
compounds (VOCs) used in these experiments were obtained either from Scott Gases or Apel- 
Riemer Environmental, Inc. The mixing ratios of the NMHC and VOC standards were between 
10 and 20 ppmv. 

2.2 Ion-Molecule Reactor 

The CIMS apparatus consists of two sections, an ion-molecule reactor where the 
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sampling and chemical ionization is performed and a vacuum system housing a quadruple mass 
spectrometer and an electron multiplier detector. Two mass spectrometers were used in this 
work. The laboratory experiments were performed on a dual channel diffusion pump system 
previously described by [Eisele and Tanner , 1993], The ambient air measurements were 
performed on the top channel of the 4 channel NCAR CIMS system as described by Mauldin et 
al. [2001]. In both cases the same ion-molecule reactor was used. 

The ion-molecule reactor consists of two sections, the transverse ion source, conceptually 
similar to that described by Eisele and Hanson [2000], and the transport tube, as shown in figure 
1. The ion source is a 7 cm x 8 cm x 2.5 cm teflon block with a 1.9 cm diameter hole machined 
through the center. A 9 cm long by 1 .9 cm inner diameter teflon tube fits snugly into the teflon 
block mating with the center hole to form the sampling inlet. Ambient air or dry N 2 is sampled 
through the inlet at 4 SLPM. A radioactive ion source, 241 Am, and the transport tube are located 
on opposite sides of the sampling tube. Between the radioactive source and transport tube is the 
ion-molecule reaction region (see Figure 1). The stainless steel transport tube is 22.86 cm long 
with a diameter of 0.762 cm and has a 0.025 cm diameter orifice located at its top. A buffer gas 
of dry N 2 is continuously flowed in front of the orifice to minimize any ambient air or un-ionized 
reagent gas from entering the transport tube. 

The two advantages of using this type of ion-molecule reactor on an aircraft platform are: 
1 ) that the transport tube provides a convenient way to extend the sampling inlet outside of an 
airplane’s boundary layer; and 2) contact between the inlet wall and the ambient air sample is 
minimized by aligning the sampling tube with the ambient air flow, thus, allowing for straight 
path sampling into the reaction region. Minimal wall interactions should also increase 
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instrument response time by inhibiting reversible desorption from the walls. 

23 Reaction Time 

The protonated ethanol cluster ions are produced by flowing 1 SLPM of an ethanol/N2 
mixture over an M1 Am radioactive ion source. 241 Am releases a particles that form (C 2 H 5 OH)H + 
through a series of reactions. These ions rapidly cluster with ethanol to form an equilibrium 
distribution of (C 2 H 5 OH) n H + . The (C 2 H 5 OH) B H + ions are accelerated across the reaction region 
through the path of the ambient sample by an electrical field. The electric field is produced by 
applying a potential of 1 .6 kV to the radioactive ion source and 100 V to the transport tube. 

The drift time (i.e., reaction time) of the ions across the reaction region can be calculated 

as 



where L is the length of the reaction region, V is the electrical potential across the region, and k 
is the mobility of the ion. However, there are two difficulties to applying equation (3) to the set- 
up used here. First, k is unknown because the ethanol cluster size distribution in the reaction 
region is impossible to predict due to a lack of thermodynamic parameters. Second, this formula 
determines the drift time for ions between two infinitely long parallel plates and does not take 
into account the specific geometries of the ion-molecule reactor or any dielectric effects from the 
teflon housing to the electric field. Therefore, an independent determination of the reaction time 
was obtained using SF 6 as the reagent ion and adding a SOj/N 2 mixture to the ion-molecule 

reactor to estimate the effective electric field strength. 

The rate constant for the SF 6 ' + S0 2 reaction has been well established as 1.0 x 10‘ 9 cm 3 
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molecs' 1 s' 1 , independent of temperature and pressure [Arnold et al. , 2000; Ferguson , 1976; Huey 
et al., 1995; Streit, 1982], The reaction time was calculated from measuring the decay of SF 6 ' as 
a function of the S0 2 concentration. From these experiments a reaction time of 7 ms was 
obtained. Using this reaction time, the reaction region length (L = 1.9 cm), and 1.90 cm 2 V s as 
the SF 6 ion mobility [Lakdawala and Moruzzi, 1980], the effective electric field strength (V eff ) is 
calculated to be 270 V. 

The largest uncertainty with estimating the reaction time lies in the calculation of the ion 
mobility of the ethanol cluster ion distribution. The relationship between ion mobility and mass 
has been difficult to quantify (e.g., [Bohringer et al., 1987; Tammet, 1995]). Along with mass, 
ion structure is known to affect mobility [Kim and Spangler, 1985]. Again, since the protonated 
ethanol cluster size distribution is unknown, it is even more difficult to relate mass and ion 
structure to mobility of the ethanol cluster ion distribution. Makela et al. [1996] measured the 
mobility of acetone cluster ions produced in a manner similar to the ethanol cluster ions in this 
study (i.e., flowing 1000 ppmv of acetone over a radioactive ion source) at 1.3 - 1.5 cm 2 V 1 s' 1 . 
Since the mass of acetone is similar to that of ethanol (58 amu to 46 amu), it is reasonable to 
assume that the ion mobility of the protonated ethanol cluster ions in this study is similar. 
Furthermore, data from Kilpatrick [1971] suggests that the ion mobility only changes 
approximately by a factor of 2 for an order of magnitude change in mass. Therefore, using 1 .4 
cm 2 V 1 s' 1 , the middle of the range from the acetone cluster study, as the ion mobility, the 
reaction time is estimated from equation (3) at 10 ms. Given the results of Kilpatrick [1971], this 
estimate is probably good within a factor of 2. Based alone on the mass difference between 
acetone and ethanol, the uncertainty in the reaction time is estimated at ± 20 - 30 %. 


June 18, 2001 


9 


23 Transport Tube Effects 

Since water from the ambient sample and ethanol from the reagent gas are prevented 
from entering the transport tube by the use of N 2 as a buffer gas in front of the transport tube 
orifice, the distribution of (CjHjOHyT cluster ions in the transport tube is dramatically altered 
from that in the reaction region. Weakly bound cluster ions dissociate in the transport tube 
allowing the (C 2 H J OH) n H + ions to reach a new equilibrium distribution. The minimum cluster 
bond strength needed to transverse the transport tube has been estimated here for a neutral, x, 
and an ion, y + , reacting to form an equilibrium cluster xy + 

X + y+ Xy+ (4) 

In this set-up, with the pressure in the transport tube at approximately 50 torr and a ~500 seem 
mass flow through the transport tube, the resulting residence time is approximately 50 ms. From 
this residence time, the largest allowable k, is 20 s '. If lq is assumed to be greater than 10 s , 
equals 5 x 10 12 . Under room temperature conditions and using 6-7 keal/mol as a reasonable 
assumption for TAS, solving for AH gives a minimum AH„ n to traverse the transport tube of 20- 
21 kcal/mol. Thus, the transport tube greatly affects what is detected since clusters with a bond 
energy of less than 20-21 kcal/mol will dissociate in the transport tube before being detected. 

As a result, the transport tube makes the observed mass spectrum less complicated and easier to 
interpret. 

Unfortunately, due to the transport tube set-up used in this work, it is impossible to know 
the detailed size distribution of ethanol clusters in the reaction region. However, it is known that 
the observed ethanol cluster distribution is quite different from that in the reaction region. 
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Because of this, in order to qualitatively test for reactivity as a function of cluster size, 
experiments were performed using ethanol/N 2 reagent gas mixtures of 24 ppbv, 900 ppmv, and 
1% C 2 H 5 OH to produce the reagent ions. In all cases, the dominant (C 2 H J OH) n H + clusters 
observed in the mass spectra were the monomer, (CjfyOItyff, the dimer, , and the 

trimer, (C 2 H J OH) 3 H + . Since the AH^ of 

(QHjOH^IT + C 2 H 5 OH -(C 2 H j OH) 3 H + (5) 
has been reported as -21.4 kcal/mol [Feng and Lifshitz, 1995], it is not surprising that no 
evidence of (C 2 H 5 OH) 4 ir or larger is seen in the mass spectra (i.e., the larger clusters dissociate 
in the transport tube before detection). Under the 24 ppbv ethanol conditions, the monomer and 
dimer were predominantly observed. At the higher ethanol levels, the monomer and dimer 
signals were seen increasing and the trimer was also observed. These results suggest that there 
are more large ethanol clusters in the reaction region under the higher ethanol conditions than 
under the 24 ppbv ethanol conditions. 

3.0 Results 

3.1 Effective Rate Coefficients 

Effective rate coefficients were determined for the ion molecule reaction of (C 2 H 5 OH) n H + 
with DMSO, NHj, CH<, DMS, and select NMHCs and VOCs using 24 ppbv, 900 ppmv, and 1% 
ethanol/N 2 reagent gas mixtures. The effective rate coefficients were determined by monitoring 
the depletion of the total ethanol signal (i.e., the sum of the monomer, dimer, and trimer signals) 
as a function of the neutral reactant added to a dry N 2 sample flow. Given the combined 
uncertainties in the estimated reaction time and the gas standard dilutions, the uncertainty in the 
effective ion-molecule rate coefficients reported here is estimated at a factor of 2. 
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The measured effective rate coefficients for the ion molecule reactions of (C 2 H 5 OH) n H + 
are reported in Tables 1 and 2 along with the proton affinities of the neutral reactants. Listed in 
Table 2 are all species tested that react with (C 2 H 5 OH) n H + at < 5 x 10' 12 cm 3 molec 1 s 1 under all 
ethanol conditions. Calculated collision rate coefficients, determined by the method of Su and 
Chesnavich [1982], in the limit that the ethanol cluster size approaches infinity, are reported for 
selected reactions. As shown in Table 1, DMSO, NH 3 , acetone, and methyl vinyl ketone (MVK) 
react with small protonated ethanol clusters (i.e., 24 ppbv ethanol/N 2 reagent gas mixture) at or 
near the calculated collision rate. The effective rate coefficients are: 1.8 x 10‘ 9 , 1.5 x 10' 9 , 1.0 x 
10’, and 1.6 x 10' 9 cm 3 molec' 1 s' 1 , respectively. DMS0 2 was also found to react with small 
protonated ethanol clusters at or near the gas kinetic rate, however, the additional uncertainty in 
the DMS0 2 standard has prevented a more precise determination to be made. Only DMSO and 
NH 3 react with the larger ethanol cluster ion distribution (i.e., 900 ppmv ethanol/N 2 reagent gas 
mixture) with effective rate coefficients of 1.5 x 10' 9 and 0.7 x 10 9 cm 3 molec' 1 s *. DMSO and 
NH 3 are also the only species to react with the largest ethanol cluster ion distribution (i.e,l% 
ethanol/N 2 reagent gas mixture), however, the DMSO effective ion molecule rate has slowed to 
0.35 x 10' 9 cm 3 molec' 1 s' 1 , while only a lower limit on the NH 3 reaction can be given because the 
system equilibrates. 

(C 2 H j OH) n H + + NH 3 - (NH 3 )(C 2 H j OH),H + + nC 2 H,OH (6) 

The ratio of product ions to reactant ions, ^ H jOH ) n h + » IS tem P erature and etha no 

concentration dependant. In the laboratory at 21 0 C and 1% ethanol/N 2 reagent gas mixture this 
ratio is 2.5 x 10' 12 . These experiments clearly show that the reactivity of the (C 2 H 5 OH) n H + ions 
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is a function of cluster size. 


3.2 Relative Humidity Tests 

The (C 2 H 5 OH) n H + ion distribution was studied as a function of relative humidity (RH). 

At low ethanol concentrations (i.e., 24 ppbv ethanol) the (C 2 H 5 OH) n H + ions are converted to the 
water clusters (H 2 0) 2 H + (37 amu) and (H 2 0) 3 H + (55 amu) for RH > 50%. But for the same RH 
under the 900 ppmv and 1% ethanol conditions, the dominant peaks are the ethanol monomer, 
dimer, and trimer. At RH < 25 % the ethanol ion peaks dominate independent of ethanol 
concentration. Thus, these experiments suggest that when the RH > 50 % it is necessary to use a 
reagent gas mixture with greater than 900 ppmv of ethanol to preserve the (C 2 H 5 OH) n H 
chemistry. It is important to maintain the (C 2 H 5 OH) n H + chemistry because (H 2 0) m H is more 
reactive than (C 2 H ; OH) D H + and shifting the ion chemistry to (H 2 0) m H + would decrease the 
selectivity of this technique. 

Additional experiments were performed to examine the effect of H 2 0 on the effective ion 
molecule rate constant for DMSO and NH 3 . This was done with a 900 ppmv ethanol/N 2 ethanol 
reagent gas mixture by adding a constant amount of DMSO, 70 pptv, or NH 3 , 400 pptv, while 
humidifying the UHP zero air sample stream from 0 to 85% relative humidity. Under these 
conditions less than 20% of the total ethanol signal was converted to water clusters, thereby, 
preserving the (C 2 H 5 OH) n H* chemistry. The effective ion molecule rate coefficients for DMSO 
and NH 3 changed by less then 20% from 0 to 80 % RH. 

3 J DMS Interferences 

The species benzene (C 6 H 6 ) and methanol (CH 3 OH) could give ion products (C 6 H 6 IT or 
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CH 3 OH C 2 H 5 OHH + ) that have the same mass as DMSOH\ However, since neither QH* or 
CH 3 OH were observed to react with (C 2 H 5 OH) n H + , they can not be considered as significant 
interferences in the detection of DMSO. Recalling that DMSO is produced by the oxidation of 
DMS, a potential interference is artifact DMSO production by ion molecule processes from 
ambient DMS present in the reaction region. To test for this possibility, DMS was added to a 
UHP zero air sample flow in concentrations ranging from 600 pptv to 69 ppbv under the 
following conditions: low humidity, high humidity, low humidity with 200 ppbv 0 3 , and high 
humidity with 200 ppbv 0 3 . The results showed that DMS did not react with (C 2 H 5 OH) n H + 
under any of the above conditions to produce a product at the mass of DMSH + , DMSOH , or C 
DMSOH + . 

3.4 Ambient Air Tests 

The protonated ethanol ion chemistry was used to measure DMSO on one channel of the 
NCAR four channel CIMS system deployed on the NASA P-3B aircraft during the GTE-PEM 
Tropics B field program in March/April of 1999. To maximize the selectivity during these 
measurements a 1 % ethanol reagent gas mixture was used. The sensitivity of the instrument to 
DMSO was measured to be 1 count per second (Hz) per pptv of DMSO at a total ethanol signal 
of 10 kHz relative to the calibration C 13 DMSO. This agrees quite well with the sensitivity of 0.9 
Hz estimated from the experimentally determined effective rate coefficient under the 1 % ethanol 
conditions. The background signal at the DMSO product ion mass was found to be no greater 
than 1 Hz (i.e., 1 pptv) which gives a limit of detection for DMSO of 0.50 pptv for a 60 s 
integration period. As expected, the DMSO sensitivity was not affected by ambient relative 
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humidity under high ethanol conditions to within ± 20 %. 

Figure 2 compares a mass spectrum taken in the boundary layer at 75% RH during flight 
16 of PEM-Tropics B, panel A, with one taken in the laboratory at 80 % RH, panel B. Note that 
the ion signal axis is linear from 0 - 1000 Hz and logarithmic above 1000 Hz. The discrepancy 
in total ion signal between the two spectra is primarily a result of the difference in pumping 
speeds between the field and laboratory instruments. As seen in figure 2, the ethanol cluster 
peaks dominate both spectra indicating the preservation of the (C 2 HjOH) n H* ion chemistry under 
these high RH conditions. Also seen in panel B are NH 4 + and (CjHjOHJFPNlV due to NH 3 
contamination in the deionized water used to adjust the relative humidity. 

Figure 3 presents a time series plot of DMSO taken during the midday boundary and 
buffer layer sampling runs during flight 16 pf PEM-Tropics B. The data shown here were 
recorded at 1 5 second time resolution. This time series highlights the response time of the 
technique. Here one can easily see the sharp change in DMSO levels when switching between 
the boundary and buffer layers. This suggests that the response time of this technique is faster 
than 1 5 seconds making it highly suitable for aircraft studies. The data also show significant 
levels of DMSO in both the boundary and buffer layers. As expected, boundary layer DMSO 
levels were 2-3 times greater than those observed in the buffer layer. A more detailed 
comparison of these observations with model simulations is presented by Nowak et al. [2001]. 

4.0 Summary and Discussion 

Overall, the protonated ethanol cluster ion chemistry presented in this work appears to be 
a highly selective and sensitive method for detecting certain atmospheric trace gases. However, 
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controlling cluster size is critical for utilizing this ion chemistry in field studies. As shown in 
Table 2, many species with PA > C 2 H 5 OH do not react with any size (C 2 H 5 OH) 0 H + clusters. As 
indicated by Table 1, for species that do react the effective ion molecule rate coefficient 
decreases with increasing cluster size. Not only does cluster size influence the effective ion 
molecule rate coefficient, cluster size also establishes the ambient conditions under which the ion 
chemistry will be effective. As discussed in Section 3.2, at high levels of water vapor small 
(Cj^OH^tT clusters are transformed into (H 2 0) m ir via 

(QHjOHyT + mH 2 0 -(H 2 0)JT + nC 2 H 5 OH. (7) 

At high relative humidity a minimum cluster size is needed to preserve the (C 2 H ; OH) d 1T ion 
chemistry and the selectivity of the technique. Thus, at high relative humidity a minimum 
cluster size is needed to preserve the (C 2 H 5 OH) n H + ion chemistry. This means that the ethanol 
cluster ion size determines the selectivity and sensitivity for a given RH range; and so, to 
effectively apply the protonated ethanol cluster ion chemistry to field measurements, the ethanol 
cluster ion size must be matched not only to the atmospheric species of interest but also to the 
ambient conditions. 

Of the species examined in this study the (C 2 H 5 OH) c H + ion chemistry is best suited for 
the detection of DMSO. Since DMSO is one of the few species that reacts with large 
(QHjOHXjfT clusters, a high level of selectivity is inherent in this ion chemistry. The use of 
large ethanol clusters also eliminates any sensitivity dependence on humidity seen in other CIMS 
techniques due to cluster ions with different numbers of water molecules reacting at different 
rates [Seeley et al, 1997]. As demonstrated under ambient conditions during PEM-Tropics B, 
the combination of high selectivity, insensitivity to ambient RH levels, and fast time response 
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and resolution makes this technique well suited for ambient measurements of DMSO. 

For many of the same reasons, this technique also shows considerable promise for the 
detection of NH 3 . However, more attention needs to be paid to the (C 2 H 5 OH)„H + cluster ion size 
than for DMSO. At the largest (C 2 H 5 OH) n H + clusters used in this study, NH 3 reaches an 
equilibrium with the ethanol cluster ions which is dependent on temperature and ethanol 
concentrations. Consequently, the sensitivity could be strongly dependant on temperature and 
ethanol concentration of the reagent gas mixture. No species in this study, expect for DMSO, 
react with (C 2 H 5 OH)„H + when using a 900 ppmv reagent gas mixture. Therefore, careful control 
of the (C 2 H,OH) n H + cluster ion size distribution is critical to maintain the high selectivity while 
eliminating the difficulties of the equilibrium reaction. In addition, sampling and background 
issues need to be further addressed as highlight by the presence of NH 3 related peaks shown in 
figure 2 panel B. 

Laboratory results also suggest that this method could be used to detect acetone and 
MVK under low ambient relative humidity conditions. Acetone is believed to be an upper 
tropospheric HO x source and to partition reactive nitrogen into PAN [Arnold et al, 1997; 
Crawford et al. , 1999; Singh et al., 2000; Singh et al, 1995]. The low water levels in the upper 
troposphere may allow for fast-time resolution acetone measurements using small protonated 
ethanol clusters. The OH-initiated oxidation of isoprene has been found to yield MVK, 
methacrolein (MACR), and formeldehyde as major products [Paulson et al., 1992; Tuazon and 
Atkinson, 1990]. Since MACR does not react with small (C 2 H 5 OH) n H* clusters (see Table 2), 
this technique could potentially be used to differentiate between MVK and MACR when 
studying isoprene oxidation mechanisms. However, since isoprene oxidation is of most interest 
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in the boundary layer, ambient relative humidity conditions are often too high for use of the 
small (C 2 H 5 OH) n H + clusters. Further laboratory study examining the possibility of using the 
small (C 2 H 5 OH) 0 H + cluster in conjunction of with a dryer to make boundary layer MVK 
measurements is needed. 

5.0 Conclusions 

The reactivity of (C 2 H 5 OH)„H + cluster ions has been found to be inversely proportional to 
cluster size, i.e., the larger clusters react with fewer species. Only two species, DMSO and NH„ 
were found to react with the large (C 2 H 5 OH) n H + cluster ions. DMSO reacts at or near the 
collisional rate with all size ethanol clusters examined. NH 3 , also, reacts with all size ethanol 
clusters, though, it reaches an equilibrium with the larger clusters. For ambient measurements at 
RH > 50%, larger (C 2 H 5 OH) n FT cluster ions must be utilized. Therefore, a high level of 
selectivity is inherent in the (C 2 H 5 OH) n H + cluster ion chemistry. Interference testing shows that 
DMSO is not produced from DMS by ion processes in the ion-molecule reactor. (C 2 H 5 OH) c H* 
cluster ion chemistry has been successfully used to measure DMSO from an aircraft platform 
during PEM-Tropics B. Laboratory studies suggest this ion chemistry be used to detect NH 3 
under ambient atmospheric conditions. 

The ion-molecule reactor described here has shown to be viable for aircraft trace gas 
measurements. While the transport tube may limit the number of species that can be detected, 
due to the dissociation of weakly bound clusters, it eliminates possible interferences from those 
same clusters, makes the mass spectrum easier to interpret, and extends the sampling inlet 
outside the airplane’s boundary layer. 
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Table l. Species found to reart with (C J H J OH) I1 H l with 
ion molecule rate coefficients > 2 x 10 12 cm 3 molec s 


Species 

P A. 


Km 

K* 


DMSO 

211.4 

2.2 

1.8 

1.5 

0.35 

NH, 

204 

1.9 

IS 

0.7 

> .0025* 

Acetone 

194 

1.9 

1.0 

nr 

nr 

MVK 

199.5 


1.6 

nr 

nr 

DMSOj 

? 


0.1 - 1.0 

nr 

nr 


PA. - proton affinity (kcal/mol) [Hunter and Lias,, 2000] 
k^, - calculated collisional rate constant 

(10 9 cm 3 molec 1 s' 1 ) [Su and Chesnavich, 1982] 

- 24 ppbv ethanol/Nj reagent mixture 
(10 9 cm 3 molec' 1 s') 

k J ia - 900 ppmv ethanol/Nj reagent mixture 
(10' 9 cm 3 molec' 1 s' 1 ) 

k„. . -1% ethanol/N, reagent mixture (10 9 cm 3 molec s ) 
nr - rate constant < 2x10 cm molec s 
* - system equilibrates 


Table 2. Species found to react with (C 2 H 5 OH) n H with 
ion molecule rate coefficients < 2x10 n cm 3 molec s 


Species 

PA. 

Species 

PA. 

Species 

PA. 

Methane 

129.9 

Trans-2-butcne 

m 

2-pen tanone 

m 

Acetylene 

153.3 

n -butane 

3 

Ethyl acetate 


Ethylene 

162.6 

lsobutane 

ns 

Cyclobutenc 


Ethane 

142.5 

Isopropanol 

189.5 

a-pinene 

D 

Methanol 

180.3 

Isobutylene 

m 

Ethyl Acetate 

D 

Propyne 

179 

Methylcyclohexane 


1 ,3-butadiene 

QR 

Propylene 

179.6 

Methacrolein 


m- xylene 


DMS 

198.6 

Cyclopentane 

? 

o-xylene 


Isoprene 

189.5 

Methylcyclopcntanc 

? 

Ethylbenzene 

i 

Toluene 

187.4 

4-methyl- 2-pentene 

? 

Methylbutenol 

D 

1 -butene 


3 -me thy 1-2 -pen tene 

194.3 

Methyl Isobutyl 
Ketone 

B 

cis-2-butene 


2-methyl-2-pentene 

ID 

II 

u 


P.A. - proton affinity (kcal/mol) [ Hunter and Lias> 2000] 
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Figure Captions 

Figure 1. Diagram of the ion-molecule reactor used in this study. 

Figure 2. Panel A, mass spectrum taken in the boundary layer at 75 % relative humidity during 
PEM-Tropics B flight 16. Panel B, mass spectrum taken in the laboratory at 80 % relative 
humidity. (Note, the y-axis is linear from 0 - 1000 Hz and logarithmic above 1000 Hz.) 
Dominant peaks in both spectrum at high relative humidity are ethanol monomer and dimer. In 
panel B, the peak present at 61 amu is isopropanol impurity present in the ethanol used to 
generate the reagent ions. 

Figure 3. DMSO time series plot for midday boundary and buffer layer sampling during PEM- 
Tropics B flight 16. DMSO observations are represented by (+) and the altitude by the solid 
line. Data are presented at a time resolution of 1 5 seconds. 
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Airborne Observations of DMSO, DMS, and OH at Marine 
Tropical Latitudes 

J.B. Nowak 1 , D. D. Davis 1 , G. Chen 1 , F. L. Eisele 1 ', R. L. Mauldin III 2 , D. J. 
Tanner 2 , C. Cantrell 2 , E. Kosciuch 2 , A. Bandy 2 , D. Thornton 3 , and A. Clarke 4 


Abstract - This paper reports the first simultaneous fast-time 
resolution measurements of dimethylsulfide (DMS), 
dimethylsulfoxide (DMSO), and the hydroxyl radical (OH) at 
manne tropical latitudes. These observations were recorded 
during the PEM-Tropics B field program on NASA’s P-3B 
aircraft. The observations of DMSO, using a Selected Ion 
Chemical Ionization Mass Spectrometry (SICIMS) technique, 
are of particular significance. They have revealed two unique 
findings: 1) average midday-tropical levels of DMSO are 
significantly higher than those predicted from current models 
when constrained by observed DMS and OH levels (e.g., 10 
pptv versus l to 3 pptv); and 2) DMSO concentration profiles 
are significantly out-of-phase with model predictions, 
maximum values being seen under near dark conditions and 
minimum values being observed at midday. Although no 
simple explanation has yet been tound tor these unusual 
results, the fact that no evidence points to a problem in the 
measurements suggests that others may exist. Clearly, if the 
observations are correct, they indicate that at least for tropical 
upwellmg regions the atmospheric sulfur budget may need to 
be adjusted to accommodate additional sources or DMSO. 

Introduction 

Since the mid-1950s, dimethylsulfide (DMS) has been 
recognized as an important species in the scientific community’s 
assessment of global climate change. This reflects the fact that 
the oxidation of ocean released DMS has been found to be an 
important marine source of sulfur dioxide (SO : ) in remote 
regions. Thus, it can ultimately lead to the formation of sulfate 
aerosol particles [e.g., Charlson et a/., 1987; Berresheim et ai . , 
1995 and references therein]. In past efforts to better quantify the 
details of the DMS oxidation process, laboratory as well as field 
studies have revealed that DMSO is a significant intermediate 
oxidation product (i.e., see Figure 1 and Turmpseed et ai .. 
[1996], Davis et ai [1998]). In fact, cumene thinking is that the 
OH/ DMS reaction is the dominant source of atmospheric DMSO. 

Here we report the first simultaneous fast-time resolution 
measurements of DMSO and the closely coupled chemical 
species DMS and OH. These airborne measurements were 
recorded over the tropical Pacific Ocean during the NASA GTE 
PEM-Tropics B field program. Because of their time resolution, 
these observations have provided an excellent opportunity to 
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explore the detailed chemistry controlling manne atmospheric 
levels of DMSO. 

Experimental Techniques and Model 
Description 

In the current study, DMSO measurements were made using a 
Selected Ion Chemical Ionization Mass Spectrometry (SICIMS) 
technique modified from that used by Berresheim et al. [1998]. 
Typically, 4 SLPM of ambient air is sampled through a v «” 
Teflon inlet tube into the reaction region where DMSO was 
converted into DMSOH" by reaction with the reagent ion 
(C:H 5 OH) x H"\ The (C : H<OH) x H* ions were generated using a 
radioactive source and accelerated across the reaction region by 
electric fields. Product DMSOH* ions were then sampled with a 
quadrupole mass spectrometer and detected using an electron 
multiplier. 

The SICIMS instrument was calibrated by the constant 
addition of a known amount of isotopically labeled C 1 ' DMSO 
( l3 CH 3 OS i3 CH 3 ) from a permeation tube. The permeation rate 
was determined gravimetncally by the manufacturer and is 
traceable to NIST standards. The tube was kept in a temperature 
controlled oven with constant N : flow maintained over it. The 
nitrogen/C 13 DMSO mixture was diluted using a muin-stage 
dynamic dilution system before being added to the inlet where 
final dilution occurred with ambient air. The typical levels ofC i3 
DMSO in the inlet were 100 pptv. Background measurements 
were performed by over filling the inlet with dry N ; or zero air. 
A potential interference considered in this system was the 
formation of DMSO from ambient DMS. Further laboratory 
testing under simulated manne conditions, however, showed no 
artifact DMSO production and no interference at the O 3 
DMSOH" calibration mass. The combined uncertainty of the 
permeation tube and dilution system is estimated at - 50% and 
primarily reflects the uncertainty associated with the emission 
rate of the DMSO permeation tube. The detection limits for 1 
and 10 minute integration times were estimated to be 1 and 0 3 
pptv, respectively, for a signal to noise ratio of 2 to l. (A more 
detailed description of the DMSO SICIMS system is provided by 
J.B. Nowak et al., Chemical Ionization Mass Spectrometry 
Technique for the Detection of Dimethyisuifoxide and Ammonia. 
submitted 2001) 



Figure 1. Abbreviated DMS/OH atmospheric oxidation scheme. 
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(pptv) 
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(pptv) 


1 Fit 15 
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- 

; Fit 1 6 

115 

7.2 -i7 9 

269 

179 - 35S 


'Fit 1" 

5.3 

2.9 - 4 4 
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114-150 



Mean 

Range 

Mean 

Range 



fpptv) 

(pprv ) 

fpptv » 

(pptv) 


1 FU 15 

3-6 

<10- 6.3 

32 

61 -96 

■Z 

‘Fit 16 

5 4 

4.2 -7.6 

108 

63 - 191 


'Fit i? 

3.7 

<1.0 -Tl 

75 

68-35 

-» -mm 


Fit 15: 16’S/U8 a W, Fit 16: 13''S/143 


Measurements of DMS were made using gas-chromatography 
mass-spectrometrv with cryogenic trapping previously described 
bv Bandy et al [1993]. The data collection rate for the PEM- 
Tropics B system was typically one sample every 5 min with a 
limit of detection (LOD) of 2 pprv. The uncertainty was =10% 
when DMS was at levels greater than five times the LOD. 
Measurements of OH were made using a SICIMS technique 
previously described by Mauldin et al. [1999]. The nominal time 
resolution for this system was 10 s although data were typically 
blocked as 1 to 5 mm averages resulting in an LOD tor this 
species in the range of 1 — 2 x 1 0 5 molec/cm 3 . The uncertainty in 
these measurements when the signal was five times greater than 
the LOD is estimated at ± 40% (2a). 

The model simulations reported here are based on a sulfur box 
model that contained 14 sulfur reactions describing both the 
sources and sinks of DMS, sulfur dioxide (SO ; ), sulfunc acid 
(H : S0 4 ), methane sulfonic acid (MSA), DMSO, dimethyl sulfone 
(DMSO*), methane sulfonate (MS), methane sulfmic acid 
(MSI A), and non-sea-salt sulfate (NSS) [Davis et ai , 1998, 
1999]. It has been used in this study to simulate the diel profiles 
of DMSO, SO;, and H : S0 4 . A critical parameter assumed in the 
current sulfur model is the conversion efficiency of DMS to 
DMSO. For the abstraction branch this efficiency was taken to 
be zero. For the addition channel, which is only 25% of the totai. 
values from 0.5 to 1.0 were used. This range reflects both the 
most recent analysis of sulfur Held data by Davis et al. [199b. 
1999] as well as DMS laboratory gas kinetic studies [Turmpseed 
ei ai. y 1996, Hynes et a/., 1993; Arsene et al y 1999]. In the 
current investigation, DMSO simulations were constrained by 
DMS and OH observations as well as by the total aerosol surface 
area. The latter was evaluated using data from a differential 
mobility analyzer (0.02< Dp <0.5pm) and an optical particle 
counter (0.148 <Dp <6pm). 

For time periods when OH data were not available, the OH 
profile was generated from a coupled photochemical time- 
dependent box model, previously described by Crawford et al 
[1999]. This model includes full HO/NO^CHa/NMHC 
chemistry. Input consists of measured values for NO, CO, H : 0, 
0 3 , UV irradiance, and nonmethane hydrocarbons (NMHCs) as 
well as for meteorological parameters [Raper et al, 2000]. 


Observations and Discussion 

DMSO was measured on flights 15, 16, and 17 during PEM- 
Tropics B. Table l presents a summary of both the DMSO and 
DMS data from these flights. Information is provided for both 
the boundary layer (BL) and the overlying transition zone 
between the BL and free troposphere (labeled here the buffer 
layer (BuL) as defined by Russell et ai [1998]). Of the BL 
DMSO observations shown, those recorded during flight 16 were 
the most comprehensive in that covered the largest fraction of the 
solar window. They also encompass the highest values recorded. 
All flight 16 BL observations were recorded at/or below 300m; 
whereas, most of the BuL measurements were made at 1500m. 


As discussed in the text below, flight 16 was the only flight 
that was downwind of the Marquesas Islands. Although the 
observations from flight 15 were geographically further 
downwind from those recorded on flight 16, they were made near 
sunset. Thus, they were some of the lowest values observed 
during the three flights in which DMSO was measured. Flight 17 
was a transit flight from Tahiti to Hawaii and the BL/BuL data 
reported were those recorded near Christmas Island at 2.5°N 
around solar high noon. Earlier flights in this region revealed 
that, on average, this location had lower midday DMS values 
than those recorded during flight 16. 

The marine region sampled dunng flight 16 has been shown to 
be strongly influenced by an island mass-effect due to the 
presence of the Marquesas Islands to the east [ Signonm el al. t 
1999]. As a result, large enhancements in surface chlorophyll 
and DMS (previously reported by Yvon et al. [1996]) are 
typically found as the Southern Equatorial Current flows around 
the islands. The net effect of the latter is that overhead 
atmospheric levels of DMS typically have been found to reach 
several hundred pptv [e.g., Yvon et al . , 1996]. 

Because of its uniqueness in providing DMS and DMSO 
levels as a function of solar intensity, the remainder of this text 
primarily focuses on the analysis of the flight 16 data As shown 
in Figure 2a BL data was recorded dunng flight 16 from 
approximately 0700 to 1330 on 7 Apnl 1999, covenng a 
significant fraction of the solar window. The sampling pattern 
executed dunng flight 16 followed an Eulenan scheme in which 



Figure 2. Diurnal trends of DMS and DMSO dunng flight 16: 
BL (bottom panel), BuL (top panel). All observational data have 
been averaged to 10 mm intervals. The la deviation iOr a i m j n 
DMSO average was = 16% tor the BL and = 25% tor the BuL. 
The la deviation in the OH data was estimated at = 19% for the 
BL and = 12% for the BuL. Similar estimates for DMS gave a 
lo deviation of ± 6% for the BL and = 12% for the BuL. The 
branching ratio and the overall lifetime tor DMSO for each 
model simulation are given in Table 2. 



Table 2. Summary of DMSO Model Simulations 


Model 

Run 

Branching 

Ratio 

DMSO 

Aerosol Lifetime ' 

DMSO 

Overall Lifetime 2 

A 

1.0 

1 hr 

0.S hrs 

B 

1.0 

6 hrs 

2.4 hrs 

C 

0.8 

6 hrs 

2.4 hrs 

D 

1.0 

2 hrs 

! 2 hrs 

E 

1.0 

13.5 hrs 

2.7 hrs 

F 

0.3 

13.5 hrs 

2.7 hrs 


Lifetime with respect to aerosol scavenging 

Lifetime with respect to aerosol scavenging and photochemical 

oxidation 


four designated altitudes were sampled multiple times but at 
approximately the same geographical coordinates (e.g., 13°S and 
143°W). The aircraft's horizontal racetrack shaped flight legs 
were flown orthogonal to the wind field and the corresponding 
fetch for this wind regime was well over 800 km. 

As seen in Figure 2, two unique aspects of flight 16 are: 1) 
levels of DMSO are quite high (e.g., BL, 7-18 pptv), and 2) the 
trend in DMSO does not follow what one might predict based on 
the DMS/OH reaction being the major source of DMSO. The 
highest observed levels are found at/or near dark conditions 
(sunrise); while, the lowest values are recorded around high 
noon. The BL DMSO trend also parallels that observed for 
DMS. As shown in figure 2a, the model simulated DMS profile, 
as constrained by a model estimated DMS flux and observed OH, 
closelv tracks the experimentally observed DMS profile. This 
agreement argues strongly against the possibility that the DMS 
and DMSO profiles were somehow controlled by abrupt shifts in 
the sampled air mass during flight 16. 

The three model profiles presented in Figure 2a, are all based 
on the assumption that the DMS'OH reaction is the sole source of 
DMSO. They suggest that the maximum value of DMSO should 
occur around 1000 LST. Model A assumes a conversion 
efficiency for DMS to DMSO of 1.0 (e.g., Arsene et at. [1999]). 
It also assumes that aerosol scavenging ot DMSO is primarily 
controlled bv its sticking coefficient (i.e, irreversible) as reported 
by De Bmyn et ai [1994]. Model B, uses the same DMS to 
DMSO conversion yield, but the scavenging etticiencv ot DMSO 
is limited by the rate of aqueous phase chemistry that forms MS 
from DMSO. In model C, it is assumed that the DMS to DMSO 
conversion efficiency is closer to that estimated from field studies 
(e.g.. 0.8) as reported in the analysis by Da\is et ai [199S] and 
that the scavenging rate is the same as used in model B Given 
this range of values for critical DMSO^ parameters, the 
scavenging lifetime was estimated to be 1 hr tor model A and 6 
hrs for both models B and C. However, upon folding in the loss 
of DMSO due to photochemistry, the overall DMSO lifetime is 
found to drop to 0.8 hrs for model A and 2 A hrs for both models 
B and C (e.g.. see Table 2). The authors note that if the 
assumption is made that the conversion efficiency for DMS to 
DMSO is closer to 0.5 (as per the laboratory kinetic results of 
Turmp seed et ai [1996] and Hynes et ai [1993]), predicted 
DMSO levels still drop another 20-30%. 

Shown m Figure 2b are the flight 16 DMSO data recorded in 
the BuL. Again, a similar diurnal trend in DMSO is seen as that 
shown m Figure 2a. The major difference between the BL and 
BuL DMSO profiles appears in the average mixing ratio values. 
For the BuL rhe range is 4 to S pptv, with an average ot :A pptv, 
whereas, for the BL it is 7 to IS pptv with an average of Ho 

pptv . 

The DMSO model simulations for the BuL are also shown in 
Figure 2b as model runs D. E. and F. The general assumptions 
staled for model runs A, B. and C, apply to the equivalent BuL 
runs D E, and F. However, because ot the .ewer aerosol oading 
of the BuL. the estimated DMSO lifetime has been increased. 
For model run D, it was 2 hrs: and for both model runs E and F it 


was 13.5 hrs. Folding in photochemical loss, the overall lifetimes 
for the above cases were 1.2 and 2.7 hrs (see Table 2). If the 
most favorable set of assumptions is made in the BuL 
simulations, the model predicted value ot DMSO can be moved 
to within a factor of 2 of the observations. This difference is 
within the combined uncertainties ot the measurements and 
model calculations. On the other hand, for the BL case, the 
difference between observations and :he most favorable model 
profile is still a factor of 3 lower at midday and over an order of 
magnitude lower for early morning hours. The largest 
uncertainty in this simulation is in the DMSO aerosol scavenging 
rate. However, even assuming no aerosol scavenging/dry 
deposition loss, the model simulated DMSO profile is still lower 
than the observations for early morning hours by a factor ot 6. 

The magnitude of the difference separating the model 
predictions and observations would initially suggest that either 
the measurements are in error or that there are problems with the 
model. Concerning possible modeling shortfalls, the most 
obvious of these is that there are other DMSO sources. In 
addition, there is the issue of vertical transport. The tact that the 
BuL had nearly half as much DMS as the BL is a strong 
indication that convection was quite significant in the -4 hours 
preceding flight 16. Both real time photographical records and a 
detailed meteorological analysis of the region using satellite 
imagerv have confirmed this. Thus, it stands to reason that the 
same convective processes that were responsible for transporting 
significant amounts of DMS from the BL to the BuL also were 
responsible for transporting a significant fraction ot BL DMSO 
into the BuL. The most effective time period for this would have 
been under dark conditions. Given this scenario, a DMS model 
simulation was earned out in which parameterized transport was 
included. The resulting DMS mixing parameter was then 
emploved for assessing the impact of transported DMSO on the 
BuL. Although this addition increased the BuL DMSO by 
another 40%, it leaves the BL DMSO discrepancy just as large as 
it was in the beginning. 

Returning to the possibility that other DMSO sources might 
exist m addition to the DMS/OH reaction, rhe tropical DMSO 
data reported by Bandy et ai [1996] prov e interesting. These are 
the onlv other tropical data having a time resolution ot less than 
an hour. Quite significant here were the findings that DMSO 
levels were also highly elevated (e.g., 1^-50 PP n ‘ Als0 - when 
all data were plotted on a common 24 hr time axis, the diurnal 
trend showed the highest values to be under dark conditions, 
decreasing toward mid-dav. Using the same suitur model, 
simulations based on measured DMS were tound to under-predict 
the observed DMSO by factors of 25 to 50. depending on the 
assumptions made concerning the DMS-to-DMS O conversion 
efficiency and the aerosol loss rate for DMSO [Chen ei ai. 
2000]. Collectively, both tropical DMSO data sets would seem 
to point toward an alternative tropical DMSO source 

Here we speculate on two other DMSO sources m the tropical 
marine troposphere The first of these involves the halogen 
reaction. BrO - DMS Previously. Tounu [1994] estimated that 
as much as 10% of the global oxidation of DMS could occur via 
the reaction with BrO. More recently, Ingham et ai. [1999] have 
reported a rate coefficient for this reaction that wcu;d double this 
prediction. In spite of these favorable kinetic developments, 
there remains two serious problems with this hypothesis. First a 
significant tropical marine source of BrO has yet to tre observe 
and second, model simulations reveal that at night BrO woui e 
rapidly depleted b> its reaction with tormaldenyce iC.H : U). 
Since, m rhe tropics, significant levels ot CH ; G persist mroug 
the mcht. the latter shortcoming would see:r.i..u > 
impossible to account for the observed daik or near aark U lb j 
observations that have been reported by liatut\ e: ... o • J an 

in the current study. cn . 

The second possibility might involve release o: DMbO from 
the ocean itself. Recognizing that a similar diurnal .<en is oun 
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between DMS and DMSO in both the Bandy data and during 
flight 16, it is tempting to speculate that these two species may 
have a common source. This possibility becomes even more 
intriguing when it is recognized that marine observations have 
shown that DMSO can be 10 to 20 times higher in concentration 
in the water column than DMS [Lee ei ai. 1999 and references 
therein]. On the other hand, since DMSO is far more soluble 
than DMS le g., Henry’s Law constant 1.17 x 10 5 versus 0.49 
mol kg' 1 atm' 1 [Watts and Brunblecombe, 1987; Sander, 2000]), 
it is difficult to see how this mechanism might explain the high 
levels of DMSO in the BL. 

Revisiting the possibility that major errors might exist in the 
observations themselves has turned up no leads. Critical here 
have been the laboratory tests which have shown no evidence of 
a DMS to DMSO conversion within the ion-molecule reactor. 
Although quite limited in number, other non-tropicai DMSO 
observations do not appear to show the same trends observed in 
the tropics [e.g., Davis et ai, 1998; Sciare et al 2000]. Both the 
presence of upwelling and deep convection at the tropical 
latitudes of this study may be important factors, but at this time 
we have not been able to show a simple causal relationship. 
Clearly, additional BL tropical field observations of DMSO, 
DMS, and OH together with water column measurements of 
DMSO are needed. If there is indeed an additional source of 
DMSO in the tropics, the natural global sulfur budget may need 
to be reexamined. 
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